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Because of the extreme complications arising from a direct deductive approach, the theory of 
the liquid state requires the use of a model involving simplifying assumptions. In this paper an 
attempt is made to formulate general principles which any such model must follow. The first step 
is a general discussion of ‘‘communal” entropy, arising from the sharing of available space by all 
the atoms. Arguments are advanced to support the contention that for a gas of hard elastic 
spheres the communal entropy is fully excited in each direction of space, and amounts in all to 
3R per mole. The communal entropy of assemblages of atoms exerting normal attractive and 
repulsive forces (in particular, the Debye solid) is considered. The geometry, the equation of 
state, and the partition function for an assemblage of hard elastic spheres are considered in 
detail. By extension of these ideas, allowing for the type of force actually exerted on each other 
by real atoms, a general form of partition function for a monatomic liquid is set up. This parti- 
tion function involves a sum of two parts, one corresponding to a vibrational motion expressed 
in terms of the Debye characteristic temperature of the solid, and the other being a translational 
term, each part carrying with it its own communal entropy. 





ECENT attempts to elucidate the prop- 
erties of liquids and the melting phenomena 

may be divided into three classes. First there is 
the attempt to proceed deductively from the 
general principles of statistical mechanics to find 
the connection between the properties of indi- 
vidual atoms or molecules and their hehavior en 
masse. This procedure is exemplified by the work 
of Kirkwood and Monroe! and that of Mayer.’ 
Another method is to set up a model, i.e., make 
certain assumptions a priori regarding the parti- 
tion function for the liquid, and use this to 
describe, rather than deduce, its properties. To 
this class belongs the work of Lennard-Jones and 
Devonshire, Alter, Walter and Eyring, Harasima, 


askin G. Kirkwood and E. Monroe, J. Chem. Phys. 9, 514 
3J. E. Mayer, J. Chem. Phys. 10, 629 (1942). 


and others. My own work‘ might be said to 
belong to this class, since in some cases models 
were used to get partition functions, but on the 
whole could perhaps be better described as 
inductive, as I have attempted to start with the 
known properties of the liquid, and from these 
infer the behavior of the molecules which give 
rise to these properties. 

Each one of these procedures has its advan- 
tages, and the papers of various authors supple- 
ment each other in such a way that we are 
gradually acquiring an understanding of the 
liquid state. That each one of these procedures 

8 J. E. Lennard-Jones and A. F. Devonshire, Proc. Roy. 
Soc. A163, 53 (1937); A165, 1 (1938); A169, 317 (1939): 
A170, 464 (1939). W. Alter, J. Chem. Phys. 5, 577 (1937). 
J. Walter and H. Eyring, J. Chem. Phys. 9, 393 (1941). 
A. Harasima, Proc. Phys. Math. Soc. Japan 21, 156 (1939). 


40. K. Rice, J. Chem. Phys. 6, 472 (1938); 7, 883 (1939); 
9, 121 (1941). 
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also has its disadvantages is obvious from the 
general agreement that, in spite of much progress, 
our understanding of the liquid state is still far 
from being satisfactory or complete. The direct 
deduction from statistical mechanics is theo- 
retically the most desirable method of procedure, 
and would answer all questions regarding the 
liquid state, if it could be carried out. Unfortu- 
nately, it leads to such great mathematical 
difficulties that it becomes necessary to make 
many approximations along the way. On account 
of the mathematical complexity, it is difficult to 
evaluate the importance of these approximations, 
and it is equally difficult to evaluate the effects 
of special assumptions regarding such quantities 
as the parameters determining interatomic forces. 
In the end, one is practically reduced to the use 
of a model. On the other hand, in the theories 
which do use a model from the beginning, there 
is often some difficulty in seeing just how this 
model is related to the actual motions and 
properties of the individual molecules. Even 
though the models may reproduce experimental 
results over at least a limited range of conditions, 
there is an atmosphere of unreality about them. 
It appears to me that what is needed is a model 
which is as pictorial as possible, and is based, as 
closely as may be, upon the actual behavior of the 
atoms or molecules in a liquid. On the other hand, 
we also require some discussion to elucidate the 
general principles which may guide us in the 
choice of a model, and help us avoid concealed 
inconsistencies. The present paper represents an 
attempt at such a discussion. In a later publica- 
tion the results will be applied to the experimental 
material available in the case of argon. 

Before proceeding to details, we indicate in a 
general way the questions we believe need to be 
resolved, and the general method of attack. The 
first matter which needs to be cleared up is the 
contribution of the “communal entropy” to the 
total entropy of solids and liquids. The term ‘‘com- 
munal entropy” was introduced by Hirschfelder, 
Stevenson, and Eyring® to describe that contri- 
bution to the entropy due to the sharing of a 
volume by many molecules, and at the same time 
was believed by them to constitute the principal 
portion of the entropy of fusion in simple liquids. 


5 J. Hirschfelder, D. Stevenson, and H. Eyring, J. Chem. 
Phys. 5, 896 (1937). 
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It was the cause of many impolite words at | 


various symposia, but recently, since the jae, 1 4 


donment of the idea that it had any direct con- | 
nection with the entropy of fusion, there has been | 


4 
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little discussion of it. It is true that the question © 


could be side-stepped by a deductive approach © 
from statistical mechanical principles, but any 
model we may introduce requires the discussion 
and evaluation of the communal entropy. The 
first portion of the paper will be devoted to a 
consideration of this matter. 

With this preliminary question out of the way, | 
we shall describe some experiments with glass 
spheres, made in the hope of elucidating some of 
the complicated geometrical problems connected 
with a random arrangement of monatomic mole- 
cules. We shall then discuss in some detail the 
statistical mechanics of a gas of hard elastic 
spheres, extending the work of Tonks,® and, 
finally, we shall proceed to a discussion of the 
general form for a partition function for the liquid 
state, basing it upon a generalization, to take 
care of intermolecular forces, of the results ob- 
tained for the gas of hard elastic spheres. Certain 
details may be left to be decided by experiment, 
and this will be done in later work. 

Throughout we use the limiting high tempera- 
ture form for all partition functions, as this is 
satisfactory for almost all liquids. 


1. THE COMMUNAL ENTROPY OF A 
CONDENSED SYSTEM |; 


In the paper in which they first used the term 
“communal entropy” Hirschfelder, Stevenson, 
and Eyring‘ gave it, in general, only a very vague 
definition, as that portion of the entropy which 
arises from the communal sharing of the volume 
available to the system by all the molecules. In 
order to make this definition precise, it is necessary 
to assign to each molecule a cell, representing its 
share of the total volume, and then calculate the 
entropy of the system assuming each molecule 
confined to its own cell ; any excess of entropy of 
the actual system in which the molecules are not 
confined to their own cells is called the communal 
entropy. In the case of a gas of point molecules, 
or any sufficiently dilute gas, the manner of 
assigning the cells is very clear; the cell for each 


6 L. Tonks, Phys. Rev. 50, 955 (1936). 
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molecule is set equal to the total volume divided 
by the number of molecules. This is exactly what 
Hirschfelder, Stevenson, and Eyring did, thereby 
obtaining a very definite value of R per mole for 
the communal entropy. The danger lies in 
applying this result to other cases. In every case 
it is necessary to consider specifically the assign- 
ment of the cells to each molecule. In some cases 
it may not be possible to make any reasonable 
assignment; then the communal entropy cannot 
be defined. It does seem to be possible, however, 
to make the assignment in the case of closely 
packed hard elastic spheres, and it is of im- 
portance in determining the partition function. 
The simple case of a one-dimensional system of 
N hard spherical atoms (or molecules—we use the 
words interchangeably) of diameter o, whose 
centers are constrained to move along a straight 
line of length /, was considered by Tonks.® For 
this system he calculated the phase integral 
explicitly; his value can be put in the form 
1\(1—0:)"/N! where 0:=No/l. If this result is 
multiplied by the momentum integral and 
divided by h, and if N! is evaluated by Stirling’s 
theorem, we get for the partition function re- 
ferred to a single atom (omitting factors which 
will be negligible, on account of the large value of 
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N, when the logarithm is taken), 
Q=(2armkT/h*)'e(1/N)(1—4;). 


(1) 


If we assign a cell of length //N to each sphere, 
and assume the sphere constrained to stay inside 
its own cell [thus giving it a free distance of 
(l/N)(1—6;) ], the partition function will be 


Q= (2amkT/h?)\(l/N)(1—61). (1’) 


The factor e in Eq. (1) arises because the length / 
is shared by all the atoms, and we see that this 
assignment of cells gives a communal entropy 
which, per mole, is equal to Rln e=R, and, in 
this case, is independent of the closeness of the 
packing. 

The assignment of cells in a two- or three- 
dimensional gas of hard elastic spheres is scarcely 
more difficult if the packing is simple cubic, and 
if the spheres are sufficiently closely packed. 
When the packing is very close, the motion of the 
atoms in any direction will always be so small as 
not to interfere appreciably with motions at right 
angles; in this case the motion may be considered 
as compounded of independent one-dimensional 
motions along the mutually perpendicular, inter- 
penetrating rows of atoms. The type of cell to be 
assigned is indicated in Fig. 1a for the two- 
dimensional case. The little square shows the 
possible region of motion of the center of the 
atom, which is the important thing. In three 
dimensions it is a little cube, whose edge will be 
a—o=v'—og=v'(1—6!) where v is the volume per 
atom, @ is the mean interatomic distance, and 
6=0*/v. In view of the independence of the mo- 
tions in the three directions of space, the partition 
function will be just the cube of that in Eq. (1) 
with v!(1—6*) substituted for (J/N)(1—61). Thus 


Q=[(2rmkT/h*)*ev'(1 — 64) }*. (2) 


This means that the communal entropy pet mole 
is equal to 3R In e=3R, one k per atom for each 
degree of freedom, since if the atoms were con- 
fined to their own cells the e would not appear 
inside the bracket. I have recently called atten- 


tion to the large value of the communal entropy,’ 


but, though it was clearly implied in Tonks’ 
work,® this was apparently not realized by 
Hirschfelder, Stevenson, and Eyring in their first 


70. K. Rice, J. Chem. Phys. 10, 653 (1942). 
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paper on the subject, nor by others (including 
myself) in the subsequent discussion.® 

We must now consider the case of a close- 
packed instead of a simple cubic arrangement. 
(We here use the term “close-packed” in the 
technical sense of describing hexagonal or, usu- 
ally, in three dimensions, cubic close-packing, as 
contrasted to the term ‘‘closely packed”’ which 
merely indicates that the atoms are close to- 
gether.) In the two-dimensional case the appro- 
priate cell assignment, which avoids any over- 
lapping of spheres, is a hexagon, as shown in 
Fig. 1b (compare Tonks’ Fig. 4). The area of the 
small hexagon which gives the free area of motion 
of the center of the atom is obviously s[(a—«)/a ? 
where s is the atomic area. In three dimensions 
the cell is a dodecahedron, the free volume per 
atom being a small dodecahedron of volume 
vl (a—a)/a * =v(1—6*)’, where v is as before the 
atomic volume and @ is given in this case by 


6=0°/a® = 03/23, (3) 


since v is equal to a*/2}. (We shall hereafter, 
unless otherwise noted, consider @ to be defined 
by Eq. (3) in the three-dimensional case. For the 
one-dimensional case, we shall continue to use 
the subscripted 6:.) We now ask whether or not 
the communal entropy is excited as before, i.e., 
whether or not we can use Eq. (2) unchanged 
except that @ is given by Eq. (3). 

We may note that there are two distinct 
questions involved here. In the one-dimensional 
case, in which the communal entropy is fully 
excited, free fluctuations of the position of atoms 
(without change in their ordinal position), are 
allowed all along the line on which they move. 
Such fluctuations, which may become very large, 
are blocked to some extent by the neighboring 
rows in the close-packed lattice, even if the 
atoms in the designated row do not move 


8 Since the communal entropy of a dilute gas is only R 
per mole, one is led to suppose that the decrease from 3R 
for hard elastic spheres at high concentration must occur 
at some intermediate stage. Presumably this will happen 
when displacements can become so great that holes of 
atomic size begin to be probable, and there is much lateral 
distortion of rows and columns. But in this situation it 
becomes difficult to assign a cell to each molecule, and it 
is not possible to set > an unambiguous definition of 
communal entropy, so that a detailed discussion of the 
change in its value is useless—fortunately, for the theory 
of liquids ideas based on a fairly closely packed set of 
spheres seem to be adequate. 
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laterally out of the row at all. On the other hand, 
if an atom does move laterally out of its row, its 
free motion parallel to this row cannot be so great 
as it was before, even if all the neighboring atoms 
remain fixed in their average positions. This | 
latter effect is presumably taken care of, at least | 
approximately, by using the value of v calculated 
for the cubic close-packed lattice and by using 
Eq. (3) in the determination of @ (equivalent to 
using the dodecahedral cell). To find out about 
the excitation of communal entropy, we must 
consider the fluctuation along a row of atoms, 
assuming the centers to be constrained to stay on 
a straight line. If the fluctuations are not too 
greatly interfered with by interlocking atoms, we 
may again consider the excitation of communal 
entropy to be complete in the three dimensions of 
space, and though the geometry of the system is 
somewhat different than in the simple cubic case, 
it would appear that we could not be greatly in 
error in assuming that the dodecahedron giving 
the free volume is effectively expanded by a 
factor of e in all its linear dimensions, again giving 
a communal entropy of 3R. 

Suppose, then, we have a row of N hard elastic | 
spheres arranged in a straight line, which can 
move freely along the line except that they can- 
not penetrate into or pass each other. With a 
length of line /, of which the portion /; is to the 
left of the edge of the atom whose fluctuations are | 
being considered, when it is in its average posi- 
tion, and the portion /, to the other side (so that 
1,+/.=l—0¢), and with the diameters of the N—1 
other atoms occupying a fraction 6; of the length 
1,+/2, a simple calculation of the standard type 
(see the Appendix) gives as the approximate 
probability of a fluctuation in position of 
amount A/ 


P=P exp [—(N—1)(Al)?/2(1—61)*%l2]. (4) 


The exponent is equal to —1 when the value of 
Al is equal to 


A=((1—-6:)/(N—1) ][2hle(N—1)]}#. (5a) 
~~ the middle atom, for which };=],.=(/—«)/2 
A=[(/—o)(1—61)/(N—1) JL(N—1)/2]!, (Sb) 


(J—o)(1—0,)/N is the free length per atom.? 
About 80 percent of the fluctuations are less than 





































9 See the end of the Appendix. 
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TABLE I, 








Fluct. ratio Entropy discrep. 











\. From Eq. (5a) or (5b), however, it appears 
that \ can increase indefinitely as N increases. 
Enormous fluctuations, therefore, become ap- 
preciably probable. 

With a cubic or hexagonal close-packed ar- 
rangement the interlocking effect of the adjacent 
atoms must necessarily be to prevent these 
enormous fluctuations. The question which must 
be answered is, therefore, how important are 
large fluctuations for the full excitation of com- 
munal entropy. That they will not be so ex- 
tremely important is clear from the discussion of 
Gurney and Mott," who pointed out that for a 
finite collection of N atoms along a line the 
communal entropy is RIn (V/N!") instead of 
R, and that this approaches R, as WN increases, 
fairly rapidly. However, we need to consider for 
different values of NV, actual values of the ex- 
pression R In (V/N!/") —R, which we shall call 
the “‘entropy discrepancy,” and, simultaneously, 
the actual probable fluctuation, which can 
best be measured as the “fluctuation ratio,” 
A: (l—o)(1-#1)/N=[(N—1)/2]§N/(N—1). 
The fluctuation ratio is the ratio between the 
probable fluctuation \ and the free length per 
atom; since this depends on N we wish to see 
whether, by assuming the atoms blocked off in 
groups of N, it can be made sufficiently small 
while the entropy discrepancy is still small. 

Before we do this, however, let us make an 
observation which will strengthen our argument. 
Consider the free space available for the motion 
of even an individual atom. Tonks has shown 
that the average distance between two next- 
nearest neighbors (i.e., between a pair of atoms 
with another atom between them) is just equal to 
what it would be if all atoms were uniformly 
spaced on the line, namely, 2//N. If, therefore, we 
fixed, in an average position, all the atoms but one, 
it would move through a distance 2(1/N)(1—6;) 


10R, W. Gurney and N. F. Mott, J. Chem. Phys. 6, 222 
(1938). 


rather than merely (//N)(1—4,), as would be the 
case if the atoms were all partitioned into equal 
lengths on the line. Thus if we considered only 
the motion of the individual atoms, assuming 
each to be free to move in the average “‘field”’ of 
its neighbors, we already have a factor 2 (instead 
of e) in the partition function, and so may be said 
in this approximation to have accounted for 
about 69 percent [i.e., (RIn 2)/R] of the com- 
munal entropy. 

Extending this same idea to a row of N atoms 
in a length /, or free length /(1—6,), we see that, 
if these V atoms are simply imbedded in a longer 
row of atoms, we should add one atomic free 
length (1/N)(1—6,) to the free length /(1—6,). 
This changes Gurney and Mott’s formula 
for the communal entropy of the row to 
Rin ((N+1)/N!/¥] and also multiplies the 
fluctuation ratio by a factor (V+1)/N. 

Allowing for these corrections, we show in 
Table I the entropy discrepancy and the fluctua- 
tion ratio for a number of different values of NV. 
For a three-dimensional case, of course, the 
entropy discrepancy must be multiplied by three 
to get the total for all three degrees of freedom, 
but even allowing for this, it appears that the 
entropy discrepancy has become quite small, say 
at N=21, before the fluctuation ratio has become 
large enough to indicate that the fluctuation will 
be appreciably affected by the interlocking 
action of the atoms in adjacent rows. Of course 
there will be some blocking of the free motion of 
a group of, say, 21 atoms, but on the other hand, 
fluctuations involving any such given group and 
an adjacent group simultaneously will be by no 
means entirely cut out, since the selection of the 
particular atom which starts the group is quite 
arbitrary. Thus it appears altogether probable 
that the assumption of complete excitation of 
communal entropy is a reasonably good approxi- 
mation, not in error by more than a few tenths 
of an entropy unit at most, even in three 
dimensions. 

This conclusion is strengthened by looking at 


the matter in another way. Consider a group of 


atoms in close-packed array, as indicated in 
Fig. 2a. If in pursuance of its fluctuations, in 
order to make the most of its possibilities of com- 
munal entropy, the atom A moves to the position 
near atom B, shown in Fig. 2b, it is true that it 
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blocks the motion of D to some extent. On the 
other hand, greater free space is left for atom E. 
Also the possibility of fluctuations of atom A 
perpendicular to the row BC is not prevented if 
it takes the position indicated in Fig. 2b. Of 
course, if it moves out of the row its freedom of 
motion in the direction parallel to that row is 
somewhat less, but, as already noted, this is taken 
account of in the definitions of v and of @ in the 
three-dimensional close-packed lattice. 

Let us now consider the situation which exists 
when the atoms are not simply hard elastic 
spheres, but exert forces on each other. In this 
case it is not easy to define a cell for each atom in 
an unambiguous way. Yet I believe that we may 
define the communal entropy in an operational 
manner which is sufficiently analogous to the way 
it is defined for hard spheres. The matter is of 
interest, for in the theory of the liquid state we 
shall wish to combine partition functions of 
different types, and we should hesitate to do this 
if the excitation of communal entropy were too 
radically different in the different partition 
functions. 

We consider first the one-dimensional case. If 
the two end atoms are held fixed in position, the 
other atoms will have equilibrium positions 
equally spaced along the line. If the vibrations 
are small enough so that the oscillations of the 
system may be assumed to be harmonic, and if an 
atom is acted upon only by its nearest neighbors, 
the partition function can be worked out com- 
pletely." Let &; be the displacement of the 7th 
atom from its equilibrium position, and let &—1 
and é;,; be the displacements of its two neighbors. 
The force acting on ith atom due to these dis- 
placements may be written 6(&:-1+£:41—28,), 
where 3 is the force constant. If now the (i—1)th 
and (4+1)th atoms are held fixed at their equi- 
librium positions, the ith atom can vibrate about 
its equilibrium position with a frequency given by 


p= a-1(b/2m)}. (6) 


If, however, all the atoms are free to vibrate, 
except the two end ones which are fixed in posi- 
tion, they execute certain normal vibrations, and 
if there are N+2 atoms (including the fixed end 
ones) there will be N such vibrations, whose fre- 


11See J. E. Mayer and M. G. Mayer, Statistical Me- 
chanics (John Wiley & Sons, New York, 1940), pp. 246-248. 
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quencies can be found by giving 7m all integral 
values from 1 to N in the following expression 
(which assumes N very large) 


v=a-"(b/m)} sin [(4/2) (n/N) ]. (7) 


It will be noted that the average of the square of 
these ‘‘perturbed’”’ frequencies is equal, as it must 


. be, to 7. 


The partition function for a single atom, 
vibrating in the field of its fixed neighbors, is 
equal to kT /hi. The motion of the single atom 
with fixed neighbors is analogous to the motion 
of a single atom through the free space equal to 
2(l/N)(1—41) (i.e., the free space from neighbor 


‘to neighbor, or twice the single atom’s quota of 


free space) in the case of hard elastic spheres. We 
might expect, therefore, that the partition func- 
tion for the whole assembly of atoms would be 
obtained by multiplying kT/hi by the factor 
e/2=1.359 (i.e., multiplying RT /2h> by the 
factor e) in order to take care of the communal 
entropy. To find the actual value of the partition 
function we must find the average value of 
In (kT /hyv) for all the frequencies given by Eq. (7) 
and take the exponential of the result. The aver- 
age of In (kT /hv) is N“ J0% In (kT /hv)dn. Using 
Eq. (7), evaluating this integral, and taking its 
exponential, we find for the partition function 
21kT /ho=1.414kT/ho. We might, therefore, say 
that in this case the communal entropy is about 
In (2X1.414)/In e=104 percent excited. 

In the three-dimensional case the problem of 
finding an exact expression, such as Eq. (7), for 
the frequencies is much more complicated.’* How- 
ever, we may use the Debye distribution of fre- 
quencies with a reasonable degree of approxima- 
tion, and in this case it is known that the partition 

2 An exact value for the partition function has been 


obtained in some special cases. See M. Born and M. Brad- 
burn, Proc. Camb. Phil. Soc. 39, 100, 104, 113 (1943). 
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function per molecule is e’kT’/hvm where vm is the 
Debye maximum frequency. Also it is known that 
vm=(5/3)4?, where > is the frequency of a single 
atom vibrating in the field of its fixed neighbors," 
so the partition function becomes (3/5)*e*kT/hi 
=1.082k7/hi. In this case, therefore, the com- 
munal entropy is only about In (2X1.082)/In e 
=77 percent excited, due to the geometrical 
factors which determine the frequency distribu- 
tion. It may be noted that ¢ will also be larger 
than it would be in a one-dimensional line of 
atoms having the same law of force, on account 
of the effects of the atoms to the side of the line 
of vibration.'* This, however, cannot be properly 
said to have anything to do with the communal 
entropy, but may have considerable to do with 
the fusion process, for, as we shall see, this 
constraining effect of the side atoms is probably 
largely, if not completely, relaxed in the liquid 
state. 


2. THE GEOMETRY OF A DISORDERED ASSEMBLY 
OF SPHERES 


Since we wish to pursue as far as possible the 
analogy between a real liquid and a gas of hard 
elastic spheres, it will be desirable to have as much 
information about the geometry of the latter as 
possible. If we suppose that at very small volumes 
the spheres are arranged in a face-centered cubic 
packing, we realize that as we allow the volume 
to increase, this regular arrangement of the 
spheres will at some time break down in favor of 
a random arrangement. The question arises, what 
will be the smallest possible volume at which a 
completely random arrangement can occur. This 
is a question which does not appear to be 
answerable by conventional geometrical methods, 
but it may readily be answered by experiment. 
As a matter of fact, experiments of this nature 
have been performed, though with a different 
object in view, by Westman and Hugill. 

Our experiments consisted of pouring water 
into a known volume filled with glass beads of 
the type used in filling fractionation columns. In 
this way the fraction of the volume actually used 
by the glass beads could be determined. Since a 
similar fraction for cubic close-packing can be 

13 See O. K. Rice, J. Am. Chem. Soc. 63, 3 (1941). 

14 See reference 13, p. 4. 


1 A. E. R. Westman and H. R. Hugill, J. Am. Ceramic 
Soc. 13, 767 (1930). 
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calculated, it becomes a simple matter to de- 
termine the increase of volume on going from 
cubic close-packing to random close-packing. 

The beads used were examined, and obvious 
off-size and defective ones discarded, to get an 
assemblage of beads which were reasonably 
spherical and fairly uniform in size. In making 
the measurements it was also necessary to con- 
sider the edge effect, due to the extra amount of 
empty space at the surface of the vessel. In order 
to estimate this effect, various sized beads and 
vessels were used, ranging from about 400 cc of 
beads of 0.44-cm diameter in an 800-cc beaker, to 
about the same volume of beads of 1.6-cm 
diameter in a 600-cc beaker. The extra volume at 
the surface was found to be equal to ~0.13 
Xradius of bead Xarea of surface. This is very 
reasonable, since it can be seen that the spheres 
at the edge tend for the most part to touch the 
surface, and have the appearance of being close- 
packed over considerable areas. It can be shown 
that, if a cubic close-packed arrangement is split 
between its densest planes, an extra volume is 
included between it and a plane surface equal to 
0.11 Xradius Xarea of surface. 

As a result of these experiments, I conclude 
that the volume of a random arrangement of an 
indefinite assemblage of spheres exceeds that of 
the cubic close-packed arrangement by 15 (+3) 
percent. Westman and Hugill obtained about 
20 percent increase. 


3. THE EQUATION OF STATE OF A GAS OF HARD 
ELASTIC SPHERES 


As already noted, the equation of state of a gas 
of hard elastic spheres has been considered in 
some detail by Tonks.® For a very dilute gas 
Tonks used the equation due to Boltzmann,!® 
which gives the pressure-volume product correct 
to the third virial coefficient, namely, 


pV/RT=1+(6/v)+§(b/2)?. (8) 


Here V is the molal volume, while v is the volume 
per molecule, and 6 is van der Waals’ excluded 


‘ volume. b= 220*, where oa is the molecular diame- 


1% TL. Boltzmann, Vorlesungen tiber Gastheorie, Part II, 
§§ 48-52. Equation (8) has been extended to the fourth 
virial coefficient by Boltzmann, K. Acad. Wet. Amst., 
Proc. Sec. Sci. 1, 398 (1899); H. Happel, Ann. d. Physik 
21, 342 (1906); and R. Majumdar, Bull. Calcutta Math. 
Soc. 21, 107 (1929). 


ter. Equation (8) has also been deduced by 
Ursell,!7 by a more general method, which 
theoretically could be extended to give the com- 
plete series. However, this involves some ex- 
ceedingly difficult calculations, and has not 
actually been carried out. Tonks’ procedure was 
to use Eq. (8) for large volumes, and the pressure 
derived from Eq. (2) for small volumes. For 
intermediate volumes he simply interpolated. 
However, using the thermodynamic relationship 
(0S/0V)r=(0p/dT)y, I showed!’ that his inter- 
polation implied a close-packed gas with com- 
munal entropy more nearly equal to 3.5R than to 
3R. I, therefore, suggested a different interpola- 
tion, which appeared just as reasonable and made 
the communal entropy come out about right. The 
interpolation does, however, cover a rather wide 
gap. It, therefore, seems desirable to try to bridge 
this gap, at least in part. Even though the gas of 
hard elastic spheres is a purely artificial problem, 
the geometrical features involved are sufficiently 
near to those of real liquids, to make any progress 
in this direction valuable. 
We proceed by using the virial of Clausius,!® 


pV =3>mu?+4r=re(r), 


where u is the molecular velocity and m the 
molecular mass, the sum being taken over all 
molecules, and where r is the distance between a 
pair and ¢(r) the force of repulsion between them, 
the double sum being taken over all pairs. For 
one mole 32mu?= RT. In the case of hard elastic 
spheres there is no force between a pair unless 
they are in collision. While they are in collision an 
average force ¢ is exerted such that, if r is the 
duration of the collision, 67 = 2uu,, where u, is the 
relative velocity along the line of centers and yu 
is the reduced mass of the colliding pair (u=m/2). 
The factor 2 occurs because the relative mo- 
mentum simply changes sign during a collision ; 
therefore, the change in momentum is twice the 
momentum itself. If n(u,)du, is the number of 
collisions occurring per second with relative 
velocities between u, and u,+du,, the virial ex- 


17 See R. H. Fowler, Statistical Mechanics (The Mac- 
millan Company, 1936), second edition, pp. 245 ff. 

18Q, K. Rice, J. Chem. Phys. 10, 653 (1942). Erratum: 
the ordinate in the figure in this note should have been 
labeled R-'0(0p/T) ». 

19 See reference 17, pp. 286 ff, or Boltzmann, reference 16, 
§§ 49, 50. 
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pression becomes 
pV=RT+ uf 2uu,on(u,)du,. (9) 
0 


In this expression the distance r at which the 
force acts is always equal to, and has been re- 
placed by, c. We now have to evaluate n(u,). For 
the first steps of this calculation we follow 
Boltzmann. Let us first consider a single pair of 
spherical molecules, and ask for the probability 
that their surfaces shall be within some very small 
distance 6 of each other. If the gas is dilute we can 
give one of the molecules any position whatsoever 
in the volume V. The chance that the other one 
will lie within a distance 6 is obviously very 
closely equal to the ratio of the volume of a shell 
of radius o and thickness 6 to the whole volume, 
i.e., 420°5/V. Since each pair is independent of 
every other in this approximation, and since there 
are N(N—1)/2~ N?/2 pairs in a gas of N spheres, 
we get 270°6N?/V for the total number of pairs 
which are within a distance 6. After a time 6/u, 
all those moving toward each other with a rela- 
tive velocity u, will have collided. The number of 
collisions of such molecules per second will 
be, therefore, n(u,)du, = 2o*u,N*f(u,)du,/ V where 
f(u,)du, is the fraction having relative velocities 
in the indicated limits and account has been taken 
of the fact that half of the pairs are moving away 
from each other. Equation (9) then becomes 


pV =RT+3 (20° N?/ vf uu, f(u,)du, 
0 


=RT[1+(b/v)], (10) 


where v is the volume per molecule, since the 
integral is twice the average relative kinetic 
energy along the line of centers, and hence is 
equal to kT. Thus we derive the first two terms 
of Eq. (8). We now depart somewhat from 
Boltzmann’s methods. We must take into ac- 
count the effect of all the other molecules on the 
probability of the collision of any pair. To do this 
we need to know how the phase integral for all 
the other molecules changes when a given pair is 
in collision. Let us now place the whole system in 
a cylinder with a weighted piston which exerts on 
the system the external pressure p, necessary to 
maintain equilibrium. In order to find the phase 
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integral of all other molecules we can assume that 
we manipulate the two colliding molecules at 
will. If these two are far apart they occupy a cer- 
tain volume and carry with them exclusion 
spheres into which the other molecules cannot 
penetrate. If they are close together these 
exclusion spheres overlap. Therefore, in order to 
keep the phase integral of the rest of the spheres 
constant, the volume of the whole gas must be 
decreased by some amount AV (the contraction is 
taken as positive), depending upon the amount 
of overlapping of the exclusion spheres. Let the 
gas now expand to the original volume, keeping 
the spheres in contact. The change in the work 
function is — pAV (the pressure, of course, in all 
these processes being practically unchanged) and 
the phase integral is accordingly multiplied by a 
factor exp (PAV/kT). This means that the proba- 
bility of collision must be multiplied by this 
factor.2° Equation (10) thus becomes 


pV=RT[1+(b/v) exp (PAV/RT)]. (11) 


If the gas is dilute we may safely assume AV to 
be the overlapped volume of the exclusion 


spheres, which turns out to be'® ;°y7? or 30, hence 
pV=RT{1+(b/v) exp [3(6/v)(pV/RT)]}. (12) 
If the exponential is expanded and the solution 
for pV/RT started, it is readily seen that this 
coincides with Eq. (8), as far as the latter goes. 
If we write pV/RT=y and b/v=x, Eq. (12) 
becomes 

y=1-+< exp (2xy). (13) 
(From the definition of 6 (Eq. 3) and 8, we find 


x =2.9640.) We can solve for y as a function of x, 


which is proportional to the concentration. How- 
ever, as x and y increase, it is found that the 
second term of the right-hand side of Eq. (13) 
shortly begins to increase very much more rapidly 
than y itself, and we soon reach a point (at about 
x=0.63 or 6=0.21) where there are no further 
solutions of Eq. (13). Since for a close-packed 
system @=1 and for the randomly packed spheres 
6 is about 0.87, it is clear that this behavior does 
not correspond to reality and can only mean that 


at great densities the value of AV cannot be as ° 


great as the overlapped volume of the exclusion 
spheres of two spherical molecules. 
20 This method of treatment amounts to use of the local 


oe See J. G. Kirkwood, J. Chem. Phys. 7, 919 
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A little consideration shows that this must, 
indeed, be the case. For, as we withdraw the two 
manipulated spheres from each other, in the 
dense gas they begin to overlap the exclusion 
spheres of other molecules before their own cease 
to overlap. We may expect no further decrease in 
effective volume available to the rest of the 
molecules after the two manipulated spheres have 
been withdrawn to something like the average 
intermolecular distance in the gas. In particular, 
when the effective free volume has become zero 
and the pressure infinite Eq. (11) demands that 
AV=0. In general we may write in place of 
Eq. (13) 


y=1+<x exp (3fxy), (14) 


where f is a fraction to be determined, equal to 
the true effective value of AV divided by its 
“ideal” value 7s70°. 

To make an estimate of f we need to consider 
the situation more closely. Let us first consider 
the case in which the gas is still fairly dilute. If 
the two manipulated spheres are far apart, they 
exclude a certain volume of space equal to the 
sum of their exclusion spheres. However, the 
actual amount of space excluded is slightly re- 
duced because these exclusion spheres occasion- 
ally overlap those of other molecules. This 
overlapping is in addition to that of the exclusion 
spheres of these other molecules among each 
other, and does not replace any of it if the gas is 
sufficiently dilute so that triple collisions are rare. 

Let us calculate the average overlap of the 
exclusion sphere of a given molecule with those of 
other molecules. In a total volume V, occupied by 
N molecules, the probability that one of the 
N-—1~WN other molecules be at a distance be- 
tween r and r+dr of the given one is 42r?Ndr/V 
=4nr'dr/v. If r<20, the overlapped volume may 
readily be shown to be® $2o*—2o*r+-r'/ 12. 
Therefore the average overlap is 


2¢ 4 1 
(4/0) f (<rot—ro'r+—ar*)ridr (15) 
e \3 12 


This is the average overlap for each one of the 
manipulated spheres. 

When the two manipulated molecules are in 
contact there is still a certain average overlap 


21 See reference 17, p. 250, or Boltzmann, reference 16, 


§ 59, 
























with the exclusion spheres of other molecules. 
However, this overlap will be less than if the 
manipulated molecules are separated, for a cer- 
tain portion, 7, of their exclusion spheres is 
already overlapped. Any other molecule is pre- 
vented from entering certain regions near one of 
the manipulated molecules, and furthermore, 
overlapping of an already. overlapped exclusion 
sphere is not fully effective in increasing the 
available volume. As a rough approximation we 
shall assume that when the two manipulated 
molecules are in contact with each other, the 
average overlap is +4 of the expression given by 
Eq. (15). Contact of the manipulated molecules, 
therefore, results in a decrease of overlap with 
others of 3% times the expression (15). 2X75 
times the expression should be subtracted from 
the ‘‘ideal”’ value of AV, or 757%, in order to get 
the effective AV or fAViaeai. The factor 2 arises 
because the overlapping with other molecules of 
both of the manipulated molecules is affected. 
The value of f is readily found to be 


f=1—(48/5mve*) 
2e 4 1 
xf (<r08—no'r+—ar* Jr, (16) 
e 3 12 


where m= 16/(2 X5) =1.6. This value of m, which 
holds only for the very dilute gas, is admittedly 
an approximation even for this case on account of 
the rough estimate of the amount of extra 
overlapping when the two manipulated atoms are 
in contact ; actually, however, the results are very 
insensitive to it, and an error in this value of m of 
twofold would not be very serious, as was found 
by actual trial. 

This treatment breaks down, of course, when 
the space becomes crowded with molecules. To 
find the average overlap in a dense gas, we shall 
consider a central molecule surrounded by m 
nearest neighbors, and we shall suppose that 
these m molecules have at their disposal a space 
equal to m molecular volumes. That is to say, we 
shall suppose them to be occupying a shell of 
inner radius o and outer radius 7; given by 


(17) 
Provided the interference of the surrounding 


molecules with each other were independent of r, 
the probability of a nearest neighbor being be- 


$m(r13—o*) = mv. 
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tween 7 and r+dr(r<r,) would be 


Anrdr/$x(r13—o*) =3rdr/(ry3—o*). (18) 


(If a molecule went beyond the distance 7;, it 
would cease to be a nearest neighbor, and another 
would in general take its place.) From this equa- 
tion one could roughly calculate the average 
overlap of a pair of nearest neighbors. It is this 
average overlap of nearest neighbors which 
should be subtracted from A V idea: in the case of a 
dense gas. For if we consider the two manipulated 
molecules we see that, if they are in what may be 
called an average position as nearest neighbors, 
we may exchange one of them with another 
“‘average’’ molecule elsewhere, and in a dense gas 
the amount of space available to the other 
molecules will still not be appreciably changed. 
In the dilute gas the overlapping of the exclusion 
spheres of other molecules with those of the 
manipulated molecules takes place independently 
of the overlapping of the exclusion spheres of 
these other molecules among themse!ves. But in 
the dense gas any motion of a molecule essentially 
becomes an exchange of position with another. If 
one of the manipulated molecules is moved from 
its average position as a nearest neighbor of the 
other manipulated molecule to a distant position, 
there is more overlapping of the exclusion spheres 
of the manipulated molecules with those of other 
molecules, but less overlapping, not only of the 
exclusion spheres of the manipulated molecules 
with each other, but also of those of the other 
molecules with each other. In a dense gas the 
possible average amount of overlapping is about 
fixed. 

However, Eq. (18) needs some modification. It 
is obvious, of course, that according to it the 
average overlapping will always be less than 
AViaeai, therefore, f must always be greater than 
zero. However, when the system approaches the 
state in which the molecules are random close- 
packed, corresponding to the assemblages de- 
scribed in Section 2 above, y must approach 
infinity. This being the case, we see by Eq. (14) 
that f must approach zero, otherwise the orders 
of infinity on the two sides of Eq. (14) would be 
different. We must suppose, therefore, that the 
average overlap approaches A Videa1, Which must 
mean that the average available free volume for 
any molecule must approach zero. The reason 
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that this happens is that the neighbors just 
outside the m nearest neighbors occupy some of 
the space inside the sphere of radius 7. In fact, in 
the random close-packed arrangement they re- 
duce it to zero.” As a better approximation we 
shall use instead of Eq. (17) 


(19) 


where c is a constant to be evaluated by setting it 
equal to mv, when v is equal to the molecular 
volume of the random close-packed state and m 
the number of nearest neighbors in that state. As 
an approximation we may set under these 
conditions 


$a(71)3—o*) =mv—c, 


y= (03/24)(12/m), (20) 


where o*/2! is the volume of the cubic close- 
packed arrangement with spheres in contact. 
This assumes that the volume per molecule varies 
inversely proportionally to the number of nearest 
neighbors. This is what actually would happen if 
the average number of nearest neighbors in a 
cubic close-packed system were reduced by re- 
moving molecules, leaving holes.”* In any event, 
the matter is not very critical. From Eqs. (20) 
and (19), then 


c=120*/2! (21) 
and from Eqs. (21) and (19) 
r= (3/4a)mv— (9/2)r—1)o°. (22) 


If we introduce a quantity ¢ such that ¢=1 if 
r,<20 and ¢{=2¢/r, if r;>20, we can, by making 
use of Eq. (18) and the expression for the 
overlapped volume just before Eq. (15), write the 
average amount of overlap in the general form 


3 ri’ 7 4aro* ar® 
—— f ( — noir \ridr 
ri—od, 3 12 





rye 4 3 
aloe 
r>—o1°L3 
1 5 
+—nif3 —-atrte (23) 
24 8 


* Actually this is an approximate view of the situation, 
for it implies that the next nearest neighbors reduce f to 
zero by holding all the m nearest neighbors against the 
central molecule. Of course, it is not true that all the 
surrounding molecules actually touch the central molecule 
in the state of equilibrium; yet f is effectively zero because 
the packing is such that the amount of free space is not 
reduced by moving an average molecule. 

* Compare O, K. Rice, J. Chem. Phys. 7, 136, 883 (1939). 
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For f we get, by dividing by ;°sa0* and subtracting 
from 1, 








r>¢3 - 1 o* 16 
ry>—o? 2 13a? 5 


9ng 1 nif8 





| (24) 


where, of course, 7; is to be determined by 
Eq. (22). However, before r; is determined, m 
must be found, and m will be a function of v. For 
the random close-packed arrangement we have 
rounded off m as 10.5. This gives by Eq. (20) an 
expansion in v of 14.3 percent over the cubic 
close-packed, which agrees well with Section 2. It 
may also be readily verified that when v is very 
large, so that o* is negligible, Eqs. (16) and (24) 
coincide. Hence we may say that the result found 
for the dilute gas is equivalent to assuming that 
there are 1.6 nearest neighbors. In between these 
conditions m should have intermediate values, 
and we shall use a system of interpolation to de- 
termine them. The essence of the situation is that 
in the dilute gas the probability of the over- 
lapping of the exclusion spheres of a given pair of 
molecules is independent of the other molecules, 
while in the concentrated gas the total amount of 
overlapping approaches a sort of saturation. The 
one state, which is equivalent to setting m=1.6, 
is dependent on the fact that in that state the 
overlap is small and f=1, while the other state, 
equivalent to setting m=10.5, results from the 
overlap being large, f=0. It seems reasonable to 
get m by linear interpolation with respect to f, 
i.e., by setting 


m=1.6+(10.5—1.6)(1—f). (25) 


We then proceed as follows. From Eq. (24) we 
get f as a function of r;/o. By Eq. (25) we get m 
as a function of f and hence of 7;/o. We then use 
Eq. (22) inversely, i.e., we find v/o* as a function 
of r/c. This, then, yields a relation between f and 
v (the latter in units o*), which is what is desired. 


_ This can be expressed as a relation between f and 


6, which is exhibited graphically as the lowest 
curve in Fig. 3. 

We may now proceed, by means of a con- 
sideration of some of the properties of Eq. (14), to 
get further evidence as to the validity of this 
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curve. We may write Eq. (14) in the form 
logio (yv—1) =logio x +0.434Xifxy. (26) 


If the left-hand side of Eq. (26) is plotted as a 
function of y, we get the heavy curve shown in 
Fig. 4. If f is known as a function of v, the right- 
hand side of Eq. (26) can be plotted as a function 
of y for any value of x =b/v, giving a straight line, 
such as those shown in Fig. 4. These straight lines 
intersect the curve either in two points, in one 
point of tangency, or in none at all. It is quite 
evident that when v is very large and x is very 
small, so that the last term of Eq. (26) amounts to 
a correction term, and y is nearly equal to 1, the 
intersection giving the lower value of y is the 
correct one. On the other hand, v approaches a 
limiting smallest value and x accordingly a 
limiting highest value, while f approaches zero. 
As this happens, it is clear that the low intersec- 
tion remains finite, but the high intersection ap- 
proaches infinity. Since y must approach infinity, 
the high intersection is to be used in this case. It, 
therefore, appears that as x increases the straight 
line represented by the right-hand side of Eq. (26) 
must eventually become tangent to the curve, 
this marking the point at which the change in our 
choice of solution of Eq. (26) occurs. This sug- 
gests that we study the conditions of tangency. 

At the point of tangency the slopes of the two 
sides of Eq. (26) must be equal; therefore 


(y—1)" = éfx. (27) 
At the same time Eq. (26) must hold. If we 
eliminate y between Eqs. (26) and (27) we obtain 


a locus of f against x (or, equivalently, f against 
6). This locus is exhibited in Fig. 3, 
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Now, if the straight line determined by the 
right-hand side of Eq. (26) for a specified value 
of x or @ is tangent to the curve of the left-hand 
side, this obviously corresponds to the largest 
possible value of f for the given value of x or @. 
Hence the true curve for f vs. @ must lie nowhere 
above the locus determined by Eqs. (26) and (27). 
On the other hand, it must touch that locus at 
one point. It is seen that the curve based on the 
linear interpolation for m barely misses the locus. 
We have, therefore, made a very slight correction 
as indicated in Fig. 3. That the correction is so 
very slight may, of course, be to some extent 
fortuitous, but there does not seem to be much 
possible leeway for the corrected curve. Using the 
latter to determine f as a function of x, we have 
obtained graphical solutions of Eq. (26) in the 
manner indicated in Fig. 4. The circled points of 
intersection in Fig. 4 are shown on curve I of 
Fig. 5, after converting x to @. This curve, we 
believe, gives the equation of state of the gas of 
hard elastic spheres with a considerable degree of 
accuracy, or at least would give it if it actually 
went over to the random close-packed arrange- 
ment at small volumes. The curve for y against 
log 6 is quite insensitive to the exact form of 
Eq. (25) up to log @=1.4 or 1.5. It is sensitive 
near the point of tangency, but in this region 
Fig. 3 appears to leave little doubt as to the 
proper choice of f. 

It is, however, most likely that if we slowly 
compressed a gas of hard elastic spheres we would 
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arrive, not at the random close-packed arrange- 
ment, but at the cubic (or, what is essentially 
equivalent, the hexagonal) close-packed arrange- 
ment. At very small volumes, the equation of 
state of the latter can be calculated from the 
partition function, Eq. (2), with Eq. (3), which 
gives” 


y=(1—64)-. (28) 


This is shown as curve II in Fig. 5, the corre- 
sponding values of f to log @=1.8 being shown as 
the broken curve in Fig. 3. In order to complete 
the equation of state we must find a means to 
join curve I at large volumes to curve II at small 
volumes; this will be done by consideration of 
the entropy. 

In the case of hard elastic spheres, where the 
accessible regions of the phase space are inde- 
pendent of temperature, f must depend only on v. 
y, therefore, is independent of temperature. We 
may thus write, since y= pV/RT, 


(8S/8V)r=(8p/8T)y =Ry/V. (29) 


y in Fig. 5, therefore, gives us R~0S/d In V 
or —R"0S/d1n 6 directly. Integration from 
log @= — 2 to very small values of log 0, along 
the curve giving the true equation of state, 
including that part which connects curve I to 
curve II, should give us the change of entropy, 
which can also be calculated independently, since 
we know the entropy at both limits, that is, the 
entropy of a perfect gas on the one hand, and the 
entropy of the cubic close-packed system, 
through Eq. (2), on the other: The change of 
entropy from a volume v2 in the perfect gas 
range, to v; in the cubic close-packed range, is 


AS=—R In [v2/(ev:)(1—68)8], (30) 


for Eq. (2) would apply to a perfect gas if the 
(1—6) factor were set equal to 1, and if Q con- 
tained a factor of e instead of e’. 

Curve I of Fig. 5 is empirically very closely 
reproduced up to the point at which it crosses 
curve II by 

y =1+2.9620+5.486?+ 21.361. (31) 
Logarithmic integration of this expression from a 
very large volume v2 (small @) to a volume 2, 


*4 Compare Tonks, reference 6, Eq. (32). 





MECHANICS OF 








LIQUIDS 13 


corresponding to the point of crossing of curve I 
and curve II (this occurs at 6=0.43) yields 
AS= —R In (v2/v;) —1.96R. 
Evaluation from Eq. (30) gives 
AS = —R In (v2/v;) —2.22R, 


which differs by 0.26R or 0.52 entropy unit. It 
is, therefore, seen that the curve connecting 
curve I to curve II, in order to give the true 
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equation of state, will have to be some such curve 
as curve III of Fig. 5, which is so drawn that the 
entropy loss on following from the left first 
curve I, then curve III, then finally approaching 
curve IT will conform to Eq. (30). This composite 
curve we propose as the most reasonable ap- 
proximation for the equation of state of the gas 
of hard elastic spheres as it is compressed to a 
cubic close-packed solid. Though it is not derived 
in a strictly deductive manner, its self-con- 
sistency from a number of different points of view 
is so satisfactory as to be strongly in its favor. At 
fairly small volumes, it deviates remarkably little 
from curve IT. 

This curve differs somewhat from the curve I 
have previously suggested,'* with less theoretical 
foundation. In the first place it follows Tonks’ 


- interpolated curve quite closely to log @=1.4 


instead of only to 1.0. However, at greater values 
of log @ it departs slightly more from Tonks’ 
interpolation, rising, as noted, slightly above 
curve II, in order to maintain the validity of 
Eq. (30). 














Curve III probably corresponds to a state of 
the system in which there are small groups of 
atoms in essentially a cubic close-packed struc- 
ture, though the system as a whole is disordered. 
On compression of the gas, these groups increase 
in size and number. They enable the system to 
avoid the high pressures of curve I, but do not 
bring it all the way down to the pressures indi- 
cated by curve II until, at the very smallest 
volumes, complete order is produced. Their 
presence makes possible transition from a dis- 
ordered to an ordered state. 


4. THE PARTITION FUNCTION FOR A 
SIMPLE LIQUID 


In the light of the above discussion of idealized 
systems, we now wish to formulate a partition 
function which can be applied with some degree 
of approximation to actual monatomic?> liquids 
of the rare-gas type. There seems to be every 
reason to believe that a considerable degree of 
disorder exists in such liquids and that there is a 
variety of configurations with the same energy. 
This would indicate that a certain portion of the 
energy of motion may be described as trans- 
lational; on the other hand, the forces between 
atoms in a real liquid are strong, and to some 
considerable extent the motion will resemble that 
in solids. We may thus say that the motion is in 
part solid-like, and in part gas-like. Of course, 
this is nothing new, and has been recognized in 
most theories of liquids. However, it is believed 
that we can lay down certain principles regarding 
the way in which these two types of motion must 
be combined, and obtain a partition function 
which can be better understood, in detail, in 
terms of these general principles than those 
previously used. 

The first question to be settled is whether the 
degrees of freedom can be divided into two 
classes, some translational and some vibrational, 
or whether it is necessary to consider each degree 
of freedom as in part translational and in part 
vibrational. 

The first idea has been considered in a some- 
what special form by Brillouin?* and Lucas,”’ and 


2 The words ‘‘atom’”’ and “molecule” are used inter- 
changeably. 
% L. Brillouin, J. de phys. et rad. 7, 153 (1936). 
27 R. Lucas, J. de phys. et rad. 8, 410 (1937). 
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can be formulated in a very attractive way. Since 
the motion of the molecules of a liquid is a 
cooperative affair, as in the case of the solid, we 
naturally assume that there are certain normal 
modes of motion. Some of these would be oscilla- 
tions of the usual type; others, however, would be 
oscillations of frequency zero, namely, transla- 
tions of a complicated kind, involving the co- 
operative motion of numbers of atoms. The 
characteristic feature of the translations is that 
they are force-free motions which include mo- 
tions which stir the atoms around, producing, 
perhaps, exchanges of position without altering 
interatomic distances, at least on the average; 
free rotations of groups of atoms would be a 
special case, and, in fact, the case visualized by 
Lucas. Attractive as this idea seems to be, I 
believe that it leads to difficulties. 

Let us consider the motion, according to the 
picture just described, when the translational 
degrees of freedom are not excited. In this case, 
the addition of all the normal modes of motion 
results in each atom vibrating in the field of its 
neighbors without changing its relative position 
with respect to its neighbors. When the trans- 
lational degrees of freedom are excited, this 
picture will not be greatly altered in the length 
of time required for the average vibration; i.e., 
the surroundings of any given atom will not 
appreciably change in this time. If the sur- 
roundings of the average atom did change ap- 
preciably in this time it would mean that the 
individual modes of motion contained an element 
of both translation and vibration, whereas we 
have postulated that these constitute separate 
modes of motion. 

Let us now consider the analogy in the case of 
hard elastic spheres in a nearly close-packed 
random arrangement. Of course, strictly speaking, 
there would be no vibrational motion in this case, 
but analogous to the vibrational degrees of 
freedom would be short range motions of the 
spheres bounded by points of collision with 
neighboring spheres, without shift of relative 
position, while longer range shifts of position and 
exchanges of neighbors would correspond to the 
translational modes of motion of the real liquid. 
There would be a rapid change in entropy with 
volume as we approached the close-packed 
random arrangement, because of the relatively 


Pade 










































1ce 


we 
nal 

ila- 

be 
sla- 

co- 
“he 
hat 
no- 
ng, 
ing 
ge; 


the 
ynal 
ase, 
tion 
f its 
tion 
ans- 
this 
eth 
i.€., 
not 
sur- 
ap- 
the 
nent 
; we 
rate 


se of 
cked 
cing, 
case, 
s of 
the 
with 
ative 
. and 
» the 
quid. 
with 


cked 








STATISTICAL MECHANICS OF LIQUIDS 15 


great change in the free path corresponding to the 
“vibrational’’ degrees of freedom, and a large 
increase of pressure as the system contracted 
toward the random close-packed arrangement. 
This would correspond to a system on curve I in 
Fig. 5, which is in itself an argument against it. 

Actually, it seems inevitable that at larger 
volumes there will be considerable mixing and 
free exchange of positions of the molecules, and 
hence contributions of both a translational and 
“vibrational” type in these individual motions. 
As we have seen, we may in general expect the 
system to follow a composite of curves I, III, and 
II, there being a tendency, as the system is 
compressed, for the atoms, at first in small 
groups, to take up the true close-packed arrange- 
ment. The reason that some atoms do this is 
that it allows more room for other atoms to 
maintain their free individual motions, which 
combine the properties of ‘‘vibration” and trans- 
lation. This means that the free path of individual 
atoms is greater than one would predict from 
“vibrational’’ modes of motion, and that the 
process of exchange of positions of neighbors 
contributes to the individual free path, and so 
takes place, or at least gets well started, in a time 
comparable to the average time of vibration. It 
is true, of course, that there is a gradual loss of 
translational motion, and that the atoms in the 
close-packed groups lose it more rapidly than the 
others. But it is questionable whether the lifetime 
of such a group is so very much greater than the 
average period of vibration. Under these circum- 
stances, it seems to me that it will be preferable 
to use what amounts to an analogy of the 
Einstein theory of the solid, in which the motion 
of each atom represents an-average of what goes 
on in the liquid. 

It will be interesting to see what this point of 
view leads to in the picture of Brillouin and 
Lucas, in which the force-free modes of motion 
are assumed to have a rotatory character. We can 
suppose the liquid to be divided up into groups 
of atoms, each group having internal vibrations 


and being free to rotate. These groups must be - 


fairly small, not over three or four atoms, in order 
for there to be enough rotational degrees of 
freedom to appreciably affect the properties of 
the liquid. If such a group can rotate freely and as 
a whole for a time greater than the average period 





of vibration, then it is legitimate to assume 
normal vibrational modes of motion and other 
normal rotatory modes of motion. If the group 
splits up, or cannot rotate freely over a period of 
this order of magnitude, then the different types 
of motion must be combined in each degree of 
freedom. It is this latter alternative which we 
believe to be the correct one. 

We may suppose that the atoms will have a 
motion similar in some respects to that of an 
atom moving in a potential field such as shown in 
Fig. 6. In the figure the potential for a harmonic 
oscillation has been split in two, and a section of 
zero potential, corresponding to translation, in- 
serted in between from point A to point B. This 
type of potential has been considered several 
times before®* and it is obvious that the partition 
function will be the sum of the partition functions 
of the harmonic oscillator and the one-dimen- 
sional “‘particle-in-a-box”’ which is free to move 
between points A and B. This would give for the 
partition function per atom T/O+ (2amkT/h?)'x, 
where 9 is the characteristic temperature of the 
oscillator and x is the distance between A and B. 

However, this expression must be modified in 
several ways. In order to see how these modifi- - 
cations are to be made, it will be best to proceed 
from the partition function for the case of hard 
elastic spheres. From the above discussion of this 
case, it appears that the partition function given 
by Eq. (2) will be valid, at least within a few 
tenths of an entropy unit, even after the volume 
has become sufficiently great so that the arrange- 
ment of the atoms appears to be quite random. 
Thus it is as if each atom had an average free 
motion in a given direction for a distance of 
vi(1— 6!) = 2-ta(1 — 6) = 2-4(a—«), where a is the 
smallest interatomic distance (distance between 
equilibrium positions of nearest neighbors) for a 
cubic close-packed lattice of molecular volume v. 
The factor e per degree of freedom takes care of 
the communal entropy. In the real liquid some of 
this free space in not truly free, i.e., some of the 
motion takes place in a force field. Let us suppose 
that the part of the motion which is completely 
free corresponds to the partition function (for 


one degree of freedom) (2mmkT/h*)!2-\e(a—ap). 


28 See, e.g., reference 5. 





Then, effectively, we may say that for distances 
smaller than dp» we are in the field of force. 

It seems reasonable to suppose that the contri- 
bution to the partition function per degree of 
freedom due to the motion in the field of force can 
be given by e'T/O(ao), where O(a) is the Debye 
characteristic temperature for a cubic close- 
packed lattice with the smallest interatomic dis- 
tance do. In this form, this part of the partition 
function carries with it its own characteristic 
contribution to the communal entropy. 

However, a little reflection will show that this 
term still needs some modification, for, in the 
liquid, the atoms are acting as individuals rather 
than as elements in a lattice. As has been pointed 
out above, the value of © in a crystal is increased 
by the presence of the atoms in rows adjacent to 
one in which the particular vibration we may be 
interested in occurs.* In fact, our previous calcu- 
lations" indicate that the force resisting displace- 
ment of a given atom from its position of 
equilibrium is, at least to a first approximation, 
just twice as great in a cubic close-packed lattice, 
and, corresponding thereto, that © is 2? times 
as great, as if only the atoms in a single row were 
effective. This result depends upon the assump- 
tion that the calculation of the frequency can 
be handled by the method of small vibrations, 
which, as will be shown in a subsequent article, 
may not be a very good assumption; still it seems 
likely that the factor 2! is a reasonably good 
approximation. 

The factor 2! is closely connected with the 
arrangement of the side atoms in fixed positions. 
In the case of the liquid the position of the side 
atoms is not fixed with respect to any pair of 
“colliding’’ atoms which we may consider. The 
side atoms may be at distances at which repulsion 
occurs, or they may be at slightly larger distances 
at which attraction occurs. Furthermore, as a 
given pair of atoms approach each other, the side 
atoms can adjust themselves with some degree 
of ease to the new situation. It, therefore, ap- 
pears that the factor 2? should not appear in the 
case of the liquid; or, if we take O(ao) as the 


* Figure 2 above may be used to illustrate this point. 
Atoms like E and D, as well as atoms like C and B, con- 
tribute the force on atom A when undergoing the displace- 
ment shown. 

29 Since we consider here the collision of a single pair of 
molecules as being more or less unaffected by their other 
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characteristic temperature of the solid lattice, it 
should be divided by 2?. But since the side atoms 
still have the same space-filling property that 
they have in the case of the hard elastic spheres, 
the partition function should be divided by a 
factor, which as a first approximation may be 
taken as the same factor which divides the free 
length in the latter case, namely 2!. Thus the 
vibrational contribution to the partition function 
becomes 


e}T / O(a) (23/24) = (2e)!T/O(ao). 


We, therefore, propose as the partition func- 
tion of a monatomic liquid, per atom, the following 
expression : 


Q=[(2e)*T/ O(a») 
+(2rmkT/h?)*2-%e(a—ao) }*, (32) 


the value of a) remaining to be determined. This 
expression, of course, gives only that portion of 
the partition function which may be referred to 
the thermal energy of the system. Of course, some 
potential energy is included in the term, thermal 
energy, namely, that part of the potential energy 
which varies with the temperature when the 
parameters, a and dp are held fixed. For the work 
function (Helmholtz free energy) of the system, 


neighbors, and since in a row of atoms, each atom is 
always under the influence of two neighbors, it might be 
supposed that © for a liquid would be less by a factor 24 
than even the © for a single row of atoms. However, if 
certain interatomic distances become smaller, an equal 
number must become correspondingly larger, if the density 
is to remain fixed. It, therefore, appears that when one 
atom approaches a neighbor, we cannot neglect the effect 
of the opposite atom, even though we neglect the effect of 
the side atoms. 
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we may write 
A=-—RT InQ+E,, (33) 


where E, is a function of a and dp and is the 
potential energy which the system would have if 
it were possible for it to approach absolute zero 
and still obey the classical statistics. In order for 
this picture to be self-consistent, E, must not 
vary with the temperature explicitly. However, 
even if the condition of constant volume is im- 
posed, which means constant a, E, may depend 
on the temperature through ap» if it varies. We 
must, in general, write 


S=-—(0A/dT)y=RInQ 
+(RT/Q)[(8Q/dT)a,a0 
+ (0Q/0d0)a,7(0a0/9T)a | 
—(0E,,/0do)a(0ao/OT)a. (34) 


There are now two ways in which we may 
regard the quantity ad». The first is to consider ao 
as an additional parameter, which may, in our 
imagination, be varied at will. Thus we may 
suppose various geometrical configurations at 
constant volume, differing in coordination num- 
ber. The more the arrangement deviates from 
cubic close-packing, the smaller the coordination 
number and the smaller the nearest-neighbor 
distance, simultaneously, must become. This will 
be reflected in a decrease in do, as the smaller the 
interatomic distance the greater the free space 
between atoms and, therefore, the greater the 
translational part of the partition function. The 
thermodynamic equilibrium will occur, at any 
given temperature and atomic volume, when —A 
is a Maximum, i.e., when 

(RT/Q)(0Q/0d0)a, 7 —(0E,/da0)a=0. (35) 


Combining this with Eq. (34) we get 
S=R In Q+(RT/Q)(0Q/0T)a.ao. (36) 


The other possible point of view is to suppose 
that do is purely a geometrical quantity, de- 
termined by the volume, the size of the atoms, 
and the forces between them. This point of view 
is, perhaps, supported by the considerations in. 
the case of the hard elastic spheres, in which no 
mention is made of apparent coordination num- 
ber, but in which it appears that full use of all the 
possible positions of the atoms must be made in 
calculation of the partition function or the phase 
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integral. Of course, the situation is more compli- 
cated in the case where the atoms exert the usual 
type of force upon each other, and we may 
question whether, at given a, the value of do is 
independent of 7. At first sight it might appear 
that this is not the case, because it would be 
supposed that, at higher temperatures where 
some of the atoms are closer together because the 
positional fluctuations are higher, it would be 
possible to find a greater range of positions over 
which the energy would be constant, and that, 
therefore, a —d» should be smaller. However, this 
would mean that the translational part of the 
energy would actually take on the character of a 
vibration, for the characteristic of a vibrational 
motion is that it has a’greater range the higher 
the temperature. The aim of the expression (32) 
is to divide up the contributions to the partition 
function into a translational and a vibrational 
part; if this is to be possible, it should be possible 
to find some a» which will not vary with the 
temperature. It results in the entropy again being 
given by Eq. (36). 

We may look at this matter in a little different 
way. The first term in the bracket in Eq. (32) is 
equivalent to another term like the second one, 
except that the equivalent free distance (corre- 
sponding to a—dp» in the second term) depends 
upon the temperature. In other words, it is like 
a collection of hard spheres, but the radius of 
these spheres depends upon the temperature. The 
effect of this dependence of the effective radius 
of the ‘“‘hard’”’ sphere on the temperature is less 
the greater the volume; hence there must be a 
term, like the second in the bracket, which has a 
lesser dependence on temperature. It is not to be 
denied that the statement that a» depends only 
on a is an assumption; it is seen, however, to be 
qualitatively reasonable, and is the simplest 
assumption that can be made in actually applying 
the theory to experimental data. 

In previously published work,‘ I have adhered 
to the first of the above propositions, namely, that 
a» is a parameter. Actually I used the coordina- 
tion number as the parameter, which is equiva- 
lent. This work ran into some difficulties in 
unpublished calculations, which I could attribute 
in part toinability to take some of the geometrical 
factors into account. This would seem to point to 
the preferability of the point of view that do is 








a function of a, determined by the geometry of 
the system. Recently, I have had some success in 
coordinating those properties of liquid argon 
which are connected with the entropy, and in 
relating them to the properties of solid argon, 
using this latter point of view. It is planned to 
publish these results in a future communication. 


APPENDIX. FLUCTUATIONS IN A ROW OF HARD 
ELASTIC SPHERES 


We consider a row of N spheres of diameter o on a line 
of length 7, and select one whose fluctuations we wish to 
investigate. Assume there are N; spheres to the left of the 
selected sphere and N2 to the right. We have 


NitNet1l=N. 


The free space to the left of the selected sphere is /;, and 
to the right is /2. We have 


Litle+o=l. 


There will be a position of maximum probability for the 
selected sphere, where the force exerted by the spheres to 
the left due to their thermal agitation is equal to the force 
exerted by spheres to the right. This is determined by 
setting these forces, which are given, respectively, by*® 


fi = NikT/(1 = 4, pli 
fe = Nok 7/G _ 6, ol. 


equal to each other. This position of maximum probability 
will not necessarily be the same as the average position; 
Tonks has shown that the average positions of the spheres 
are regularly spaced along the line. The position of 


30See Tonks, reference 6, Eq. (15). Note that our 
.1= No /1; and A2= Noo /Is. When fi =fo, then 4,1=61,2=4. 
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maximum probability will, however, not be displaced 
from the average position by more than the greater of 
(L/N)(1—61:)/2N1 or (1/N)(1—0:)/2Ne. If the selected 
sphere be given a displacement of Al (assumed small com- 
pared to /; or J/2) from the position defined by setting 
fi: =fe=fo (we shall hereafter designate this position by the 
subscript zero) we have 

a NRT =f (1- — a...) 

(10h oF Al hi o(1 = 01) 7” 

NokT wee Al 

ft" G6 a i+) 

fe—fi=fo(l—o)Al/(1—41)/i, of», 0. 


The amount of work necessary to effect the displacement 
will be 
»Al 
[, (fa—fuddAl = foll—o) (Al)?/2(1 — 8,)hs, os, 0. 
The probability of a displacement Al will, therefore, be 
given by 
P =Py exp[—fo(l—o) (Al)2/2(1— Oli, oa, kT I, 


where Pp is the probability of a displacement of zero. 
Now 


— NRT _ NT 
* A=—Adle (1—-O)lee 
(Nit+N2)kT (N—1)kT 





7 (i — 61) (Li, o+/2, 0) ~ (1—6:)(l—0) 
Hence 
P=P, exp[— (N- 1 )(Al)?/2(1 — 6;)71, oe, ol], 


which is Eq. (4). It should be noted that @:, in these 
equations is equal to (V—1)@/(/1,0+/2,0) or (N—1)0/(l—@) 
rather than No/l, which means that 1—@, is equal to 
(l— No) /(l—c) instead of (I— No)/l. The free length per 
atom, for the N atoms, which is (-— No)/N, is, therefore, 
equa! to (1—61)(l—o) /N. 
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An analysis is made of the collection and definition of a beam of ions by an ion source. It is 
found that the first slit of the source has no discriminating effect but the second slit does have. 
This is due to the fact that ions of certain values of initial energy, angle of velocity vector, and 
position of striking the plane of the first slit cannot get through the second one. An efficiency of 
collection of a source for ions of a certain initial energy is defined as the ratio of the number with 
that initial energy that get through the second slit to the number with that initial energy which 
get through the first slit. An expression for this efficiency is derived and it is applied to rather 
typical data to show how an abundance ratio obtained by taking the ratio of peak heights may 
be in error several percent. The method of measuring isotopic abundances is recommended 
wherein the focusing conditions are changed only by varying the magnetic analyzer field. 


OR some time it has been believed that it is 
possible for the accelerating slit system in a 
mass spectrometer ion source to discriminate be- 
tween ions of different mass, i.e., to collect ions 
of one mass more effectively than those of 
another. This would lead to different values of 
the isotopic abundances measured with different 
instruments. 

One possible source of discrimination is the 
fact that the solid angle of collection, i.e., the 
angle in which ions are collected, may vary from 
one ion to the next during a determination of 
isotopic abundances. It is the purpose of this 
paper to analyze the reason for this discrimina- 
tion and to calculate the magnitude of its effect. 
The effect of this change of solid angle of collec- 
tion on peak shape for ions from dissociating 
molecules has been considered by Hagstrum and 
Tate. 
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1H. D. Hagstrum and J. T. Tate, Phys. Rev. 59, 354 
(1941). 


Let us consider an ion source such as shown 
schematically in Fig. 1 where the ions are formed 
by a beam of electrons coming into the ionization 
region parallel to the plane of the first slit. In 
this figure are shown two ions originating in the 
region of the electron beam with different values 
of 6, one ion being able to traverse the second 
slit and one not being able to. Let us further 
restrict ourselves to an ion source which is 
magnetic field free.* Since all ions are formed 
with an initial velocity (either thermal or due to 
dissociation energy, as for ions from molecules), 
let us denote this velocity by c and denote the 
angles of its vector by 6 and ®. Let us denote by 


/ 


c’, 6’, ®’ the corresponding quantities for the 
ion as it passes through the first slit, and by c’”’, 
6’’, &’’, the corresponding quantities as it passes 
through the second slit. Let V; be the potential 
difference between the first slit and the position 
traversed by the electron beam. Let V2 be the 
potential difference between the first and second 
slits. Let the separation between the first and 
second slits be z and the slit width in each case be 
2d. We shall consider an ion source in which 
no magnetic field is present, and we shall assume 
that the electric field is uniform both in the 
region above the first slit and in the region be- 
tween the first and second slits. 

If the area of the first slit is small compared to 
the area of the electron beam and if the first 
slit is located well away from any of the edges 
of the electron beam, it is clear that there can 





2.N. D. Coggeshall and E. B. Jordan, Rev. Sci. Inst. 14, 
125 (1943). 
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be no selective influence due to the first slit 
alone. For, consider ions of any mass, there will 
be created in the electron beam a definite number 
per unit area; and since the electric field is uni- 
form, the number arriving per unit area on the 
plane of the first slit will be equal to this, so the 
number collected by the first slit is proportional 
only to its area and the number created per unit 
area in the electron beam. 

The second slit, however, does have a selective 
influence as can be seen from Fig. 1. This is be- 
cause only part of the ions going through the 
first slit get through the second slit and this 
fraction depends upon the initial velocity of the 
ions and upon the accelerating voltages used. 

Let us consider ions of mass m. We shall de- 
note by p(c) the number of ions of mass m and 
initial velocity c which are created per second, 
per unit area in the electron beam. Now since 
the first slit gathers as many, per unit area, as 
are created per unit area in the electron beam 
we have: 

f(0’, ®’, c’)dw’ = p(c)(dw/4r), (1) 
where f(6’, &’, c’) is the number of ions of co- 
ordinates 6’, &’, and velocity c’, contained within 
unit solid angle w’ at the plane of the first slit; 
dw refers to solid angle in the electron beam. 
We assume that the ions have initial velocities 
randomly distributed. Using db=d®’ and ¢ sin 6 
=c’ sin 6’, f(6’, ®’, c’) turns out to be: 

p(c) c’? cos 6’ 


f(0, ®, c’) =— . (2) 


2e c(c?—c’? sin? 6’)} 





We may now evaluate the selective influence of 
the second slit by defining an efficiency of collec- 
tion p(c). We define p(c) as being the ratio of 
ions of mass m and initial velocity c getting 
through the second slit to the number getting 
through the first slit. From previous considera- 
tion we know that this is equal also to the ratio 
of the number of ions of mass m and initial 
velocity c in the final ion beam to the number 
created in the electron beam in an area equal to 
the area of the first slit. 


poa’—f [70 ®’, c’)dw'da’ 


p(c)= » (3) 
, p(c)A’ 


where A’ is the area of the first slit. The integral 
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is integrated over those portions of solid angle 
w’, such that ions striking the first slit in those 
portions of solid angle do not get through the 
second slit; also it is integrated over the area of 
the slit. 

In order to evaluate the above integral, let us 
consider the limiting values of 6’, 6’ that an ion 
of initial velocity c may have at the first slit and 
still get through the second slit. Let an ion 
strike the first slit, at a distance s from one edge 
and with a given ®’ and initial velocity c. Then 
one limiting value of 6’ must satisfy the following 
relationships: ; 

c’ sin & cos ®’z 
o™ ’ 
Tn, 
where /Z,, is the average value of the component 
of the ion’s velocity perpendicular to the planes 
of the slits, while the ion is between the first 
and second slits, if V2> Vj as is usually the case, 
we may write Hy =c’’/2 since 2eV2/m>c?, where 
¢ is initial thermal velocity or even dissociation 
velocity of a few volts. 
This gives: 
s=(2c’/c’’) sin # cos ®’z. (4) 





The other limiting value of 6’ is given by: 
0’ =sin-! (c/c’), (5) 


which is the case when all the initial velocity c 
is parallel to the plane of the first slit. The limits 
of ®’ are: 


sc’ 


+cos~! — 
22¢ 


so p(c) becomes, neglecting end effects: 
p(c)2Ld — f(c) 








p(c) = : (6) 
p(c)2Ld 
where L is the length of the slits and f(c) is given 
by: 
we. c’ cos—!(se’’ /2zc) 
fO= [a sf de’ 


cos 6 sin 6’d6’ 


sin—!(e/c’) 
xf 2 ot \h 
sin—!(se’’/2zc’ cos &’) (c?—c 2 sin? @ yj 





The factor 2 appearing in f(c) is the alterna- 
tive to setting up an additional integral where 
instead of having ds we have d(2d—s), i.e., 
the integral to take care of ions having ®’ differ- 
ing by 180° from the values considered. 
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Fic. 2. Energy in electron volts. 


Integrating with respect to 6’ and inserting 
the limits we obtain: 


2p(c)L 2d cos—!(se’’/2zc) 
fo-—— ff 
rg 0 Yo 
any? 3 ; 
x(1-(—) ) de'ds. 
2zc/ cos? ®’ 
Suppose we let (sc’’/2zc) =a 


f= ff 
1 




















x (cos? &’ —a*)id&'ds. 
cos &’ 
Suppose we next let 
cos? &’ —a? 
1—cos? &’ 
we obtain: 
2p(c)L r™ r® v2(1—a*)duds 
fo=—— ff 
T 0 0 (1 +v?) (v?+a?) 


which can be integrated by partial fractions to 
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yield: 
2d 
flo)=e(OL J (1—a)ds; (7) 
0 
or using value of a and inserting in (7) we get: 
dic” 
fle) =2e(e)La( 1 ~), (8) 
z2c 


Now this expression for f(c) has been derived for 
values of c such that the second slit was opera- 
tive in cutting out ions regardless of where they 
struck the first slit. However, for small values of 
c there are regions of the first slit where the ions 
may strike and be free of the discriminating 
effects of the second slit. For such values of c 
the integral f(c) becomes 


22¢ 





? 


c 


flc)=o(\L f(1—a)ds where §= 


or 


f(c) =p(c)L2(c/c’’). (9) 
The value of c for which f(c) changes from form 
(8) to form (9) is of course given by d=2c/c’’, 
i.e., it will be given by the root of the equation: 


- Qe i 
a(e+— Vit v:)) = $C. 
m 


Let us call this transition value @. Then our 
efficiency of collection is given by: 


p(c) =1—[(1/2)(2/d)(c/c’’)] for (10) 
b(c) =[(1/2)(d/2) (cc) } for (11) 


Or if we express the initial velocities in terms 
of electron volts of energy V we obtain: 


c<é, 
crc. 


z V , _ 
p(v)=1-—(——__) for V<V, (12) 
2d\V+V,4+V2 
d V+Vi4tV2\! m 
o(v)=—(——) for V>V. (13) 
22 V 


As they should, the two expressions for p( 1) 
agree in value and slope for V)=V. Graphs of 
p(V) for various values of z/d and for V;+ V2 
= 1000 volts may be seen in Fig. 2. These curves 
agree in general over-all characteristics with 
those obtained by Hagstrum and Tate.! How- 
ever, in the important region of small initial 
velocities the curves differ radically. 

As the function p(V) is independent of the 
mass of the ion being collected, true values of 
isotopic abundance ratios should be obtained if 
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all the ions are collected with voltage conditions 
remaining the same in the ion source. The only 
way this can be realized, of course, is by bringing 
first one ion beam and then another into focus 
by varying the magnetic field only. 
Let us assume that the ions when formed from 
a monatomic gas will have a Maxwellian energy 
distribution. Then the density of ions of isotope 
M,, in the path of the electron beam, as a func- 
tion of energy will be 
Pal V)dV = KpgAe** /%T Vd V, (14) 


where Pa=partial pressure of isotope WJ 


e>* 1 
A= rv(—) — , 
tkT/ (300)! 
I"=energy in electron volts, N=total density 
of atoms of all isotopes, e=electronic charge, 
m,=mass of isotope Ma. 

K represents the fraction of atoms present in 
electron beam which are ionized; i.e., it repre- 
sents ionization efficiency. We shall assume K 
to be the same for all isotopes of same element. 

The intensity of the ion current for isotope MM. 
will then be proportional to J, where: 





l= f Pal V)P(V)dV. 
0 


Since e~**/*%T<1 for V of the order 1 volt 
and for temperatures of the order of 500°C we 
may use the expression for P(V) which was ob- 
tained for small values of V, i.e., (12): 


1.=KpoA [ oe 300kT 174 
0 


Z V ; 
x(:-|—_] av. 
2d L. Vit Vot+ V. 


Since for operating conditions V«(Vi+ V2) 
we have 


1.=Kp.A f e-#¥ 30087 Vi(1 —g, VIVA, 


where 
1 1 Zz 


2+ Vs) Nd 


i! (300)*rk*7T?\ 3 
jeter 


lo 3 
300RkT *) 
“| 
é 


Za 


i (15) 





COGGESHALL 








_{l ze | 
Ig=paA |; —e ’ 
2300kT © | 


_ 300k7\? 
A= Ka(— -) ' 
é 


Suppose 7°=450°K, then 


Ta= PaA{4.46—ga}. (16) 


where 


By a similar calculation we can show that the 
ion beam intensity of an isotope MM, would be 
proportional to J, where: 


I= prA {4.46 —g5}. 


Using the usual practice of taking ratios of 
peak heights to get isotopic abundances, the 
observed abundance ratio of isotope M, to M, 
would be: 


Ta/Ty= (Pa) (4.46 — ga) / (bo) (4-46 — go) ] 


instead of p/p» which would be the true abund- 
ance ratio. 

As a test of the order of magnitude of the error 
introduced by such discrimination let us consider 
the neon isotopes Ne”? and Ne”. Let us suppose 
in our instrument that the ions are focused for a 
path with radius of curvature 5’. Further, let us 
suppose that we have a magnetic field of 800 
gauss and the two slits in the ion source are }” 
apart and 0.007” in width. It would be necessary 
to accelerate the Ne” ions through a voltage 
V,t+V2=965 volts and the Ne” ions through 
V,+ V2=877 volts. Using these values in our 
expression for J,/J, we obtain: 


I(Ne*®) p(Ne?®) 
= x 
I(Ne**) p(Ne??) 


1.02. 





In other words, the value obtained by taking a 
ratio of the peak heights would be 2 percent in 
error. This is much larger than the usual experi- 
mental error and it is certainly outside the limits 
of accuracy some investigators claim. 


DISCUSSION 


In Fig. 2 are to be seen graphs for P(V) for 
values of V up to 10 volts and for Vi+ V2=1000 
volts with different values of z/d. As can be 
seen from them, the smaller the value of 2/d 
the less discrimination an ion source will have; 
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that is, the less the separation between slits as 
compared to slit width. When working with ions 
originating in dissociation processes we would 
expect the characteristics of the ion source 
strongly to affect the shape of the peaks. The 
dependence is through the P(V) curve and is 
discussed by Hagstrum and Tate.! From our 
Fig. 2 we can see that the P(V) curve would 
strongly affect the peak shape due to the steep- 
ness of the curve for energy values of the order 
of one electron volt. 

It is to be remembered that the analysis above 
is for a rather simplified ion source, i.e., one in 
which the electric fields are uniform and one 
which is free of a magnetic field. The latter condi- 
tion is not hard to meet but the first condition is 
one which in practice could only be approxi- 
mated. This is because of the fringing of the 
fields near the slits and the penetration of the 
strong field into the weaker field region through 
the first slit. 

In view of the simplified ion source analyzed, 
the numerical results and the expression for 
P(V) should not be taken as exact for any ion 
source. However, it is reasonable to believe that 
the qualitative implications are true and that the 
expressions obtained are good as a first approxi- 
mation for any actual case. Actually it is ex- 
pected that the discrimination would be worse 
since the penetration of the high voltage field 
through the first slit can give the ions an added 





amount of velocity in a direction parallel to the 
planes of the slits and in this way decrease their 
chances of getting through the second slit. 
However, the method recommended above for 
getting abundance ratios, i.e., focusing by vary- 
ing the magnetic field only, would be expected 
to give the best answer as this would not vary 
the discrimination for the different ions as their 
peaks were measured. 

From Eq. (16) one would expect the peak 
height of any ion to grow if it were focused using 
progressively higher values for 1;+1>. This 
is illustrated in Fig. 3. In this figure not the 
peak height is plotted, but a quantity to which 
the peak height is proportional, namely, I/pA. 
In the plotting directly of the peak the slope 


‘would be greater or less depending upon units 


used. As may be seen the curve, which is for 
any ion focused for values of 1,+ 1: between 
600 and 900 volts, is slightly concave downwards 
but could be represented quite accurately by a 
straight line. 

It would be desirable to have a similar analysis 
for the type of ion source in which there is a 
magnetic field present in the ion source, usually 


. for the purpose of aligning the ionizing electrons. 


Some discussion has been made of this type of 
source® but the amount of discrimination is not 
calculated. 


*E. B. Jordan and N. D. Coggeshall, J. App. Phys. 13, 
539 (1942), 
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Viscosity and Structure of Pure Hydrocarbons 
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The two constants, A and E, of a two-constant equation, g=A exp (— E/kT), for the fluidity 
gy, are tabulated for a number of hydrocarbons. The variations of these constants with structure 


are discussed. 


INTRODUCTION 


NUMBER of investigators, among whom 
may be mentioned Thorpe and Rodger,! 
Mabery and Mathews,’ Wilson and Allibone,* 
Spilker,4 Hugel,®5 and Mikeska,* have investi- 
gated the effect on viscosity of various structural 
relationships. However, all these workers at- 
tached too much importance to the viscosity at 
some arbitrary temperature, such as room tem- 
perature, the boiling point, or the melting point. 
There has as yet been no attempt to eliminate 
completely the temperature variable in such 
investigations. 
Kauzmann and Eyring’ have shown that the 
viscosity » of a liquid can be expressed by the 
equation: 


Ash 
AoAsA? 


e ASF /ReAH*F/R A 





y= 


where ); is the distance between moving layers 
of molecules; A» is the dimension of the unit of 
flow, i.e., the molecule, in the direction of flow; 
\3 is the dimension of the unit of flow perpendicu- 
lar to \; and dg; A is the distance moved by the 
unit of flow in passing over a potential energy 
barrier; AH™ is the activation energy per mole, 
i.e., the energy required for a molecule to reach 
the top of a potential energy barrier; AS” is 
the corresponding entropy of activation; T is the 
temperature; and R is the gas constant. If we 


‘T. E. Thorpe and J. F. Rodger, Phil. Trans. A185, 397 
(1894). 

2C. F. Mabery and J. H. Mathews, J. Am. Chem. Soc. 
30, 992 (1908). 

3W. J. Wilson and B. C. Allibone, J. Inst. Petroleum 
Tech. 11, 180 (1925). 

4A. J. Spilker, Zeits. f. angew. Chemie 39, 997 (1926). 

5G. Hugel, Chemie et Ind. 26, 1282 (1932). 

6L. A. Mikeska, Ind. Eng. Chem. 28, 970 (1936). 

7W. Kauzmann and H. Eyring, J. Am. Chem. Soc. 62, 
3113 (1940). 
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write the equation in terms of fluidity ¢, since 


— 


we have 
NoA3A? 


Aih 


erS*/Re—AHF/RT. 





© 


Since (A2A3A?/A1) ~ (V/N), the molecular volume, 
the equation may be simplified to: 
V 
o = —e4S*/Re-SH*/RT, 


Nh 


Since V and AS* do not vary greatly with tem- 
perature, this equation is of the form suggested 
by Arrhenius :’ 


b= A-#IRT * 


If we consider A and E constants (experi- 
mental results justify this assumption, at least 
to within the experimental error of most vis- 
cosity measurements), these two quantities can 
be used to compare viscosities of various liquids. 
Values of E and A have been calculated for 44 
hydrocarbons by the method of least squares using 
viscosity data reported in the literature. The 
results are tabulated in Table I. 


VALIDITY OF THE ARRHENIUS EQUATION 


If A and E are constant, a straight line should 
be obtained when log ¢ is plotted against 1/T. 
Actually there is a very slight curvature, but in 
no case does logio¢ vary from the calculated 
value by more than 0.027, corresponding to a 
6.5 percent error in ¢. Differences of more than 
0.015 in log @ (3.5 percent error in 0) are rare 


8S. Arrhenius, Medd. Vetensk. Nobel. 3 [20] (1916). 
* Recently R. M. Barrer, Trans. Faraday Soc. 39, 48 
(1943) has considered this type of equation. 
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VISCOSITY OF HYDROCARBONS 
TABLE I. 

No. Name and formula A X1073 EX10-3 Reference 
£. n-Pentane, CsHi2 5.520 1.504 a 
Z n-Hexane, CeHi, 5.598 1.690 a 
3. n-Heptane, C;His 5.842 1.858 a 
4, n-Octane, CsHis 6.641 2.088 a 
5. n-Nonane, Co9Heo 8.485 2.394 b 
6. n-Decane, CioH2» 8.936 2.565 b 
7. n-Undecane, Cy,H 9.940 2.778 b 
8. n-Dodecane, Cj2Ho. 11.51 3.010 b 
9. n-Tetradecane, Ci4H30 11.40 3.203 ¢ 

10. n-Hexadecane, CigHa, 16.12 3.652 c 

11. n-Octodecane, CigH3¢ 18.40 3.891 c 

12. n-Dotriacontane, C.H¢¢ 18.66 4.922 Cc 

13. Benzene, CeHe 9.803 2.428 a 

14. Toluene, CsH;—CH; 5.777 2.059 a 

15. Ethyl-benzene, CsH;CH2CH:; 5.336 2.084 a 

16. n-Butyl-benzene, CeH;CyHg 3.139 2.033 d 

i. 2-Methyl-butane, C—~C-—t— 6.099 1.549 a 

c 
18. 2-Methyl-2-butene, C—C : C—C 4.811 1.375 a 
! 
Cc 
19. 2-Methyl-1,3-butadiene (Isoprene), C : i < 5.382 1.434 a 
C 

20. 1,5-Hexadiene, C : C—C—C—C :C 5.823 1.628 a 

71. 2-Methyl-pentane, sie ileal 5.818 1.676 a 

ee 

72. 3-Methyl-pentane, iat ‘Niall 5.194 1.638 e 

c 
C 
| 

23. 2,2-Dimethyl-butane, C—C—-C—C 8.002 1.978 

| 
& 

24. 2,3-Dimethyl-butane, a 8.563 2.035 e 

c.< 

25. 2,3-Dimethyl-1,3-butadiene, C : + << 7.054 1.724 e 

| 
S &: 
26. 2-Methyl-hexane, C—C—C—C—C—C 9.122 2.070 a 
c 
| 
27. 4,4-Dimethyl-1-pentene, C : C—C—C—C 9.057 2.194 f 
( 
28. 2,5-Dimethyl-hexane, C—C—C—C—C—C 6.142 1.974 e 
| | 
ik Cc 
29. 2.382 e 


2,7-Dimethyl-octane, i iat ii 7.113 


C 
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TABLE I.—Continued, 

















No. Name and formula A X107% E X10-3 Reference 
=e 
30. Cyclopentane, | Cc 5.554 1.878 e 
i 
c—C 
C—C 
‘ 
31. Methyl-cyclopentane, cc 6.603 2.055 e 
C—C 
c—C 
/ \ 
32. Cyclohexane, C Cc 15.21 2.905 e and a 
cc 
C—C 
33. Methyl-cyclohexane, C C—C 2.352 1.678 f and d 
‘or of 
cc 
34. 1,2-Dimethyl-cyclohexane, C ‘Cc 17.98 3.074 f 
com of 
| | 
C—C 
ct 
3S. 1,3-Dimethyl-cyclohexane, "si ‘c—C 8.150 2.353 f 
ee 
é 
c—c 
37. 1,4-Dimethyl-cyclohexane, C—C” Cet 14.89 2.861 f 
C—C’ 
38. 1,2-dimethyl-benzene, o-CgH,(CHs): 6.141 2.277 a 
C—C 
C y % 
39. 1-cyclohexyl-2-methyl-propane, C—C—C—C cc 8.787 2.422 d 
be of 
40. Triphenyl-methane, (CsH;);CH 15.38 4.600 d 
41. Biphenyl, CeHs—CeHs 32.07 4.226 d 
Cc C 
, i 
Cc c 
42. Decahydro-naphthalene, | | | 7.252 3.034 d 
i i C 
VANS 
G c 
CH CH: 
oe le. TY ale, “Oe 
CH C CH, 
43. 1,2,3,4-tetrahydro-naphthalene, | | 9.800 3.112 d 
CH Cc CH: 
ee 
CH CHe 
44, Naphthalene, CioHs 580.3 5.977 d 








* Hydrogen atoms are omitted. C—C indicates a single bond, C : C 
a double bond. 

¢ Landolt-Bérnstein, Tabellen, main volume, Part I, pp. 130, 149-52. 
>’ Bingham and Fornwalt, J. Rheol. 1, 372 (1936). 





¢ Landolt-Bérnstein, Tabellen, third supplement, p. 162. 

4 Landolt-Bérnstein, Tabellen, second supplement, pp. 100-105. 
¢ International Critical Tables, Vol. VII, pp. 216-18. 

4 Landolt-Bérnstein, Tabellen, first supplement, p. 134. 
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VISCOSITY 


and may safely be attributed to errors in a 


particular viscosity measurement. The usual 
difference is about 0.005 (1 percent). It should 
be noted that, while viscosities can be measured 
very accurately, the temperature is often not 
reported accurately, particularly in high tem- 
perature measurements, and this gives rise to 
large deviations from calculated values. In view 
of these considerations it may safely be said that 
the Arrhenius equation holds to a fair degree of 
accuracy, and the deviations may be attributed 
to variation of \;, A, AS*, and AH* with tem- 
perature. 


EFFECT OF VARIOUS STRUCTURAL FEATURES 


By comparing groups of compounds which 
differ in only one structural feature, we can de- 
termine how this feature affects A and E. Since 
A is very sensitive to errors in E, comparison of 
the A values cannot be put on a quantitative 
basis, but qualitative trends can be observed. 
Since the E values are more accurate, greater 
confidence can be placed in our conclusions 
concerning them. 

The results may be summarized as follows: 

1. Addition of a methylene group to a molecule 
increases E by about 200 calories and A by 
about 1000 poise. am 

2. A and E are both increased greatly by 
branching the carbon chain. Since branching 
causes a smaller increase in E for an equal change 
of A than does lengthening the chain, a highly 
branched hydrocarbon of low molecular weight 
should have a smaller temperature coefficient 
than a straight-chain compound of sufficiently 
high molecular weight to have the same A value. 

3. Branching seems to have a greater effect 
on compounds which are already branched than 
on normal hydrocarbons. 

4. The position of a side chain has very little 
effect on A and E. More data will be necessary 
before we will be able to say which arrangements 
tend to have higher values for these constants. 

5. A isa linear function of chain length, while 


AE becomes a little smaller with increasing . 


chain length, as is shown in Fig. 1. A redetermi- 


OF HYDROCARBONS 
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Fic. 1. Variation of A and E with chain length. 


nation of the viscosity for m—C39.H¢4 would be of 
great significance. 

6. Unsaturation seems to decrease A and E, 
but more data are necessary before the effect 
of various kinds of unsaturated bonds can be 
determined quantitatively. 

7. Cyclic compounds have much larger values 
of A and £ than the corresponding aliphatic 
compounds. In the formation of a five-membered 
ring A is increased by about 500 poise! and E 
by about 380 calories. In the formation of six- 
membered rings, A is increased by about 10,000 
poise and E by about one kilocalorie. 

8. Reduction of an aromatic ring to an ali- 
cyclic increases A by about 9000 poise! and E 
by about 700 calories. 

9. A and E seem to be very much larger for 
condensed ring systems than for single rings, but 
the data on these compounds are not very 
consistent nor complete. 

We wish to thank Professor Eyring for his 
helpful discussions. 
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A Direct Experimental Determination of the Electron Affinity of Oxygen* 
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The electron affinity of oxvgen atoms for one electron, the negative of the energy change at 


absolute zero of the reaction 


c+ Ogas—>O ~ gas» 


has been determined to be 70.8+2.0 kcal./mole. The method used was to determine the ratio of 
negative ions to electrons leaving a hot filament surface exposed to oxygen gas at a low pressure. 





GENERAL 


T is possible to determine experimentally the 

electron affinity of an atom, i.e., the negative 

of the energy change at absolute zero accom- 
panying the reaction 


, | —_ . 


provided that the number of gas atoms and the 
ratio of electrons and ions leaving a surface of 
known temperature can be determined, pro- 
vided that the accommodation coefficient is 
unity, and provided that information is avail- 
able for the calculation of the standard free 
energy of each species at unit pressure and the 
experimental temperature. This method has al- 
ready been successfully employed in the determi- 
nation of the electron affinity of chlorine? and 
iodine.? The method is extended here to the 
investigation of the electron affinity of oxygen 
for a single electron, corresponding to the reac- 
tion at absolute zero 


Qens HO get O ene (1) 


Measurements were made using oxygen gas as 
a carrier over a range of temperatures from 
2030°K to 2230°K, of pressures of oxygen from 
0.1 micron to 2 microns of mercury, and over a 
20-fold range in electron current. An average 
value of 70.8+2.0 kcal./mole was obtained, 
which value evidenced—within the limits of 


*Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

t Submitted in partial fulfillment of the requirements for 
the degree of Doctor of Philosophy at Columbia University. 

1 Kenneth J. McCallum and Joseph E. Mayer, J. Chem. 
Phys. 11, 56 (1943). 

2 J. J. Mitchell and Joseph E. Mayer, J. Chem. Phys. 8, 
282 (1940). 

’P. P. Sutton and Joseph E. Mayer, J. Chem. Phys. 3, 


20 (1935). 





experimental error—no trend with temperature, 
pressure, or electron current. 

The method of calculation of the electron 
affinity of oxygen is similar to that used by 
McCallum and Mayer for chlorine! with the 
exception that the assumption of complete 
dissociation of the molecule cannot be made in 
this case, the percent dissociation ranging from 
29 to 90. Instead, one obtains the expression 


I; My; 3 
—AEo®= RT, |n (~) 
I.Po Me 


2xm.kT.\' 2kT, Qio 
CYL a 
h? a Qio- 


where —AE,’ is the electron affinity of oxygen; 
I;/I, is the experimental ratio of ion and electron 
currents; m; and m, are the masses of the oxygen 
ion and electron, respectively ; T; is the tempera- 
ture of the filament; @ is the value of an at- 
mosphere pressure expressed in dynes/sq. cm; 
Qio and Qio- are the internal partition func- 
tions of the oxygen atom and ion, respectively ; 
and R, k, and h have their usual significance. 
The quantity Po is the pressure of a hypothetical 
oxygen atom gas at the temperature of the fila- 
ment for which the number of atoms hitting the 
filament would be equal to the number leaving 
under the conditions of the experiment. This 
pressure is to be calculated from the observed 
pressure of oxygen gas surrounding the filament 
and the dissociation constant of oxygen. 

The exact numerical value of the internal 
partition function for the oxygen ion is not 
known with certainty since its spectrum has 
not been observed experimentally. However, an 
approximate value may be obtained from the 
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following considerations. Since the electron con- 
figuration of the oxygen ion is similar to that of 
the fluorine atom, the ground state is ?P, the 
low lying term being the *P3/2 with the next 
higher, ?P 1/2, separated from it by less than the 
separation of these two terms in the fluorine 
atom or 407.0 cm~. Moreover, extrapolation of 
the curve obtained by plotting the log of the 
separation in wave numbers of the ?P3,;2. and 
the ?P1;2. terms as a function of the log of Z—2, 
where Z is the nuclear charge of the ion, of the 
isoelectronic sequence AIlV, MglIV, NalIl, 
Ne II, and FI indicates that the separation of 
these terms for the O- ion is less than 190 cm 
(see Fig. 1). Higher terms do not contribute at 
the experimental temperatures. Qio lies, there- 
fore, between 4+2e-274/7s and 6, and it can be 
assumed that 


Qio- =6.0, (3) 


without causing an appreciable error in the 
calculation, since the temperature 7, is about 
2000°K. 

The ground level of the neutral oxygen atom 
is an inverted *P, consisting of the *P2, *P;, and 
3Po terms. There are no other low lying levels 
with the higher 'D, and 'So terms lying high 
enough to be neglected. The internal partition 
function of the neutral oxygen atom is © 


Qio = 5+3e-226-2/ Ts 4 g—324.5/ 7s, 


The log of this quantity, between 2000°K and 
2400°K, may be approximated by 


log Qio =0.9124+0.0934210-‘7,. (4) 


Substituting the above expressions for Qio- 
and Qio in (2), the resulting equation for —AEo®, 
in kcal./mole, may be written 


— AEo®= (4.574T,/1000) (log A; 104/A, 
—log Po+9.963 X10-‘*T, 
—1.12110-"T,2+4.793), (5) 


in which A;/A, is the ratio of the galvanometer 
deflections due to ions and electrons, respec- 
tively; Po is the pressure of oxygen atoms in 
microns, and 7, is the experimental filament 
temperature. 

The calculation of log Po from the observed 
pressure of oxygen molecules may be accom- 
plished by the following method. The sum of the 


number of oxygen molecules leaving the filament 
plus one-half the number of oxygen atoms leav- 
ing the filament is equal to the number of oxygen 
molecules striking. The number Z of particles 
of mass m striking a unit area of surface in unit 
time is related to the pressure and temperature 
by the equation Z=P/(2rmkT)'. Substituting 
this expression for each species in the material 
balance above yields the following equation for 
the pressure of a hypothetical atomic oxygen gas 
at the temperature of the filament, from which 
as many atoms would hit the filament as are 
leaving it under the conditions of this experi- 
ment: 


Po=2-'!Kp{[1+8P/Kp(T./T,)*]'—1}. (6) 


In this Kp is the dissociation constant for the 
reaction O2:=20, T, is the temperature of the 
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Fic. 1. Plot of log Ay2p3—2p; vs. log (Z—2). 


molecular oxygen gas striking the filament (room 
temperature), 7, is the temperature of the gases 
leaving the filament (filament temperature), P 
is the observed pressure of O2 gas in atmospheres, 
and Po is the pressure of atoms in atmospheres. 

The assumption is made in this that the ac- 
commodation coefficient of the oxygen molecules 
hitting the surface is unity, i.e., that the relative 
numbers of atoms and molecules leaving the 
surface is the same as if the filament were sur- 
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rounded by a gas at filament temperature. 
In support of this assumption we may point 
out that experiments of Bryce* have shown that 
the dissociation of Hz on a tungsten filament 
follows the simple equilibrium hypothesis. This 
has been shown by calculations of Doty® based 
on Bryce’s data. 

It is interesting to note at this juncture that 
an investigation of the history of the oxygen 
atoms after leaving the hot filament indicated 
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Fic. 2. Schematic diagram of apparatus. 


that they did not persist as atoms and neither 
reacted with nor were absorbed by the metal 
parts or walls, but instead recombined to form 
oxygen molecules. The percent recombination 
was found to be 100+10. 

Returning to the calculation of Po, the dis- 
sociation constant Kp can be calculated in the 
usual fashion from statistical mechanics, em- 
ploying the necessary spectroscopic data of 
atomic and molecular oxygen. The resulting 
expression for Ko contains several exponentials, 
and can be greatly simplified by taking the 
log of this expression and expanding by Taylor 
and McLaurin series about the temperature 
T)>=2200°K. The resulting equation for log Kp, 
which may be used from 1980°K to 2420°K, is 


log Kp = —25,617/T+5.8720+4.4450 
x 10-*7—0.5669X10-772. (7) 


EXPERIMENTAL 


The apparatus used in the determination of the electron 
affinity of oxygen is shown diagrammatically in Fig. 2. 


*G. Bryce, Proc. Camb. Phil. Soc. 32, 648 (1936). 
° Paul M. Doty, J. Chem. Phys. 11, 557 (1943). 
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Electrolytic tank oxygen of better than 99.5 percent purity, 
obtained from the Ohio Chemical and Manufacturing 
Company, was slowly introduced at A and then purified by 
passing it through a liquid nitrogen trap B as a precau- 
tionary measure to remove water and the halogens, the only 
substances likely to cause an appreciable error in the 
determination of the electron affinity of oxygen. The 
mercury in the cut-off D was raised and the pressure in the 
2-liter ballast bulb and behind the capillary leak C was 
allowed to rise to such a value as to give the desired pressure 
in the experimental and control tubes, after which the 
oxygen supply was cut off by closing stopcock A. In this 
way a constant pressure was maintained behind the leak, 
the amount of gas passing through the capillary having a 
negligible effect on the pressure in the ballast bulb. After 
passing through the experimental and control tubes, which 
were in parallel, the oxygen flowed through a Princeton 
model diffusion pump with suitable backing. With this 
arrangement of introduction and pumping of gas, it was 
possible to provide a constant flow through the tubes 
(amounting, under operating conditions, to about 1300 
cc/sec.) at any pressure between 0.001 and 100 microns of 
mercury. 

The construction of the experimental tube and the ex- 
perimental procedure used for the measurement of the ion 
and electron currents were essentially the same as those 
described in a previous paper.! The tube contained a 
cylindrical plate, between similar guard rings, and a 
radially shaped grid, all being coaxial with a tungsten 
filament. Ions and electrons leaving the filament were 
accelerated to the plate by a positive potential, whereas the 
electrons could be deflected to the grid by a magnetic field 
induced by a solenoid coaxial with the filament. 

A more accurate means for obtaining the voltage across 
the filament than that described in reference 1 was used, the 
detailed construction of which is shown in Fig. 3. The upper 
and lower filament supports were welded to tungsten 
through glass leads and were constructed of General 
Electric tungsten wire with a diameter that was (18)? times 
that of the tungsten filament. (With this ratio of diameters, 
within the precision of the temperature measurements, the 
voltage correction for the cooling effect of the leads could be 
simplified by assuming that the temperature at the ends of 
the filament was one-quarter of the maximum tempera- 
ture.®) The lower support consisted of a bifilial coaxial 
spring and was provided in order to keep the filament taut 
and centered within the radially shaped grid. 

The filament current passed through one lead of each 
of the bottom and top supports while the remaining leads of 
the two supports were used to measure the potential across 
the filament. The filament itself was a 10-cm General 
Electric tungsten wire with a diameter of 0.02 cm and was 
spot-welded to the upper and lower supports. The spot 
welding caused little distortion of either the filament or 
supports at the points of contact. With this arrangement, 
the voltage across the filament could be measured directly. 
In addition, the error in the voltage due to thermal po- 


6H. Jones and I. Langmuir, Gen. Elec. Rev. 30, 310 


(1927). 
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tentials was found to be less than 0.2 percent at 2000°K, as 
determined by reversing the current through the filament 
and again measuring the potential drop. This latter 
procedure was carried out in the measurement of all tem- 
peratures and was particularly significant in the measure- 
ment of the cold resistance of the filament. 

In order to examine carefully the temperatures of the 
filament in a vacuum in the experimental tube as calculated 
from the five commonly used functions of current and 
voltage,® a control tube was put in the gas flow line in 
parallel with the experimental tube. It consisted simply of 
a filament, taken from the same piece of stock wire from 
which the filament in the experimental tube was taken, 
with similar supports and in a tube of approximately the 
same size. The filament of this tube was connected in series 
with that in the experimental tube so that the current in 
each filament would be the same. During a run in vacuum, 
the voltage across each filament and the current were 
measured, and these data were combined with the indi- 
vidual filament dimensions in order to calculate and com- 
pare the temperatures of each. It was found that for a given 
filament the temperatures as calculated from all five 
functions agreed to within 10 to 20°C in both the experi- 
mental and control tubes, with the disagreement of the 
temperatures of both filaments as calculated from the same 
function being even less than 10°C. 

The resistance function was used to calculate the filament 
temperature when runs were made in oxygen, since other 
functions give erroneous values in the presence of gas. 
Nevertheless, it is to be noted that calculation of the 
temperature by means of the other functions gave values 
which did not differ from the resistance temperature by 
more than thirty degrees, with the agreement in nearly all 
the cases being better than 20°C. Calculation of the temper- 
ature difference to be expected under the experimental 
conditions used indicates that one should not observe a 
difference greater than 10°C. 

The cause for the discrepancy found in temperatures 
calculated from the various functions was due primarily to 
the fact that it was necessary to determine the diameter of 
a filament from the measurement of the cold resistance 
before and after a run, a method which is less accurate than 
the usual method of weighing and measuring the length of 
the filament to calculate the diameter. It was necessary to 
employ this less accurate method since the filaments were 
attacked by oxygen, about five percent of the molecules 
striking a given filament reacting to form tungsten oxides 
and thus decreasing the diameter. For this reason, it was 
also necessary to limit the time of a run to about twenty 
minutes. 

The ion currents encountered in this work gave galva- 
nometer deflections of about 3 cm which were constant to 


about 5 percent. This precision was attained by surrounding, 


the experimental tube by a Dewar of liquid nitrogen, which 
reduced the fluctuations in both the ion and electron 
currents. This was presumably due to the removal of 
condensible gases given off by the filament and grid. There 
was no other apparent source of condensible gases, in- 
cluding mercury, since the tubes were protected from these 





by liquid nitrogen traps on the oxygen inlet and exit tubes 
and on the tube to the McLeod gauge used in measuring 
oxygen pressures. 

With these low ion currents it was important to eliminate 
the possibility that part of the ion current might possibly be 
due to electrons, resulting from their incomplete removal 
by the grid. The completeness of removal of electrons by 
the grid was checked in two ways. First it was found that 
varying the current in the solenoid from about 4 to 7 
amperes had no effect on the ion current, the ion current 
remaining constant. In addition, the oxygen supply was 
quickly cut off after a reading by raising the mercury in 
cut-off E (Fig. 2), and in every case, the ion current dropped 
to zero. 


DISCUSSION OF RESULTS 


In Table I are tabulated the results of sixteen 
experimental determinations of the electron 
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Fic. 3. Filament supports. 


affinity of oxygen as calculated from the indi- 
vidual sets of data, employing Eqs. (5), (6), and 
(7). The percent error for each value is also indi- 
cated and was calculated from the estimated 
precision of the individual measurements of 
temperature, pressure, electron current, and the 
ion current in each set of data. The error in 
temperature was from 2-4 percent; in pressure, 
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TABLE I. 
Estimated 
total 
O+e =O- Max. O2+e =O2- 
Kp —AE—° o —AEo° 

No. T.°K Pymicrons O2 % dissoc. O2=20 kcal./mole error Weight kcal./mole 

1 2150 0.067 81 4.71X10-° 75.62 9.5 1 83.27 

2 2230 0.37 85 1.351075 72.18 6.0 1 81.11 

3 2090 0.50 53 2.00 x 10-6 72.63 6.0 1 74.46 

4 2230 2.0 59 1.24 10-° 68.92 5.2 Z 72.04 

5 2230 2.0 59 1.2410~ 68.37 5.1 2 70.49 

6 2230 1.6 63 1.24 10-5 65.81 5.5 2 69.67 

7 2090 1.10 40 1.97 10-8 70.33 a | 3 70.01 

8 2110 0.86 48 2.61 X 107 70.51 2.8 3 71.57 

9 2220 0.50 79 1.10 10-5 69.76 3.9 3 77.11 
10 2220 0.50 79 1.10 10-5 71.05 3.6 3 78.40 
11 2180 1.15 58 6.62 X 10-6 69.88 3.8 3 72.72 
12 2180 1.07 59 6.62 X 10-* 69.82 3.9 3 72.87 
13 2030 1.00 29 8.34X 1077 73.59 S| 3 71.41 
14 2170 0.094 90 5.10 10° 66.47 8.9 1 77.07 
15 2170 0.129 87 5.10 10-* 64.66 9.9 1 74.34 
16 2170 0.129 87 5.10 10-6 65.43 9.2 1 75.10 

Average value 70.8 73.8 














10-50 percent; and in ion-electron current ratio, 
12-44 percent; resulting in a 1-3 percent, a 
0.5-4 percent, and a 1-4 percent maximum error, 
respectively, in the calculated value of the elec- 
tron affinity. The average value of the electron 
affinity was calculated by weighting each value 
inversely to the square of the percent error in 
the customary fashion. The average value 
calculated in this way was 70.8 kcal./mole with 
a standard deviation of 0.8 kcal./mole. 

The authors do not believe that the experi- 
mental certainty is as great as the formally 
calculated error of 0.8 kcal./mole would indicate. 
The tacit assumption of a zero reflection coeffi- 
cient may not be strictly true and cannot be 
corrected for since there are no experimental data 
that indicate the actual value of the reflection 
coefficient of oxygen and similar molecules. The 
assumption of a reflection coefficient greater 
than zero would increase the calculated electron 
affinity. In view of these considerations, it is 
not unreasonable to assume a larger probable 
error of 2.0 kcal./mole in the average value of 
—AEo°. 

Since the calculated election affinity should be 
independent of the conditions, such as tempera- 
ture, pressure, percent dissociation of oxygen 
molecules, or electron emission, in the reaction 
vessel, a plot of the values of —AEo° obtained 
in Table I as a function of these variables should 





result in a horizontal line. In Fig. 4 are shown the 
results of some of these plots. Similar plots of 
—AEo° as a function of the calculated value of 
the equilibrium pressure of oxygen atoms or the 
electron emission yielded curves similar to that 
against Pos. These curves are only approximate 
since they were obtained by the method of least 
squares with an arbitrary assignment of weights 
of 1, 2, or 3 to the individual points, a more 
precise weighting being unwarranted. The curves 
indicate that there is slight, if any, trend of the 
obtained electron affinity with temperature, 
pressure, electron emission, or percent dis- 
sociation. 

Several values of the electron affinity of oxy- 
gen have been obtained by other methods. The 
empirical extension of ionization potentials of 
the appropriate isoelectronic series yielded a 
value of 88 kcal./mole.’ Similar calculations have 
also yielded the values 23,* 50,9 and 36° kcal. / 
mole. A direct experimental determination of 
the electron affinity of oxygen has been made by 
Lozier and by Hanson by the analysis of nega- 
tive ion currents produced by collision of low 
energy electrons with oxygen,'' carbon mon- 

7G. Glockler, Phys. Rev. 46, 111 (1934). 

8 R. Bacher and S. Goudsmit, Phys. Rev. 59, 948 (1934). 
— and M. Maltbie, Zeits. f. Physik 75, 748 
( 10 Hellmann and M. Mamotenko, Acta Physicochimica 


7, 141 (1937). 
11 W. Lozier, Phys. Rev. 46, 268 (1934). 
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oxide," and nitric oxide.” The values of 51+5 
and 46+9 kcal./mole, respectively, were ob- 
tained. 

The reason for the disagreement between the 
values obtained by the method of Lozier and 
that obtained by this method is not apparent. 
There exists the possibility that the ions leaving 
the hot filament, under the experimental condi- 
tions used here, consisted of oxygen molecule 
ions instead of atomic oxygen ions. However, 
this seems unlikely since the electron affinity of 
molecular oxygen has been determined experi- 
mentally by two different methods and has 
yielded the values 3 kcal./mole'* and 8 kcal./ 
mole.!* On the other hand, the electron affinity 
of molecular oxygen has been estimated from 
cycles involving the heat of dissociation and 
heat of solution of HO», leading to a value of 
62 kcal./mole. 

It is to be noted at this point that the possi- 
bility of the formation of O= is excluded since 
the total electron affinity of oxygen correspond- 
ing to the reaction O+2e=O>= has been calcu- 
lated to be negative by 225 kcal./mole.'® 

Assuming that the ions observed were O.~ 
resulting from the reaction 


O.+e=O0--, 


calculation of the corresponding electron affinity 
for each set of data in Table I, yielded the values 
given in the final column of Table I. The equa- 
tion used in calculating these values was derived 
by the same method used to obtain Eq. (5), re- 
sulting in the following expression for the elec- 
tron affinity of molecular oxygen: 


—AEo®=4.574T,/1000[log.(7;x 10*/ I.) 
—log Po» (observed) +7.896 X10-‘T, 
~0.897 X10-7T,2+-4.704]. 


It was assumed in this derivation that the parti- 
tion functions of O2 and O.~ have approximately 
the same value, which is probably not in error 
by more than 2 kcal./mole. 


2G. Hanson, Phys. Rev. 48, 476 (1935). 


13N. E. Bradbury, J. Chem. Phys. 2, 827 (1934); ibid. 


2, 835 (1934). 

141. B. Loeb, Phys. Rev. 48, 684 (1935). 

18 J. Weiss, Trans. Faraday Soc. 31, 966 (1935). 

1 F, R. Bichowsky and F. Rossini, The Thermochemistry 
of the Chemical Substences (Reinhold Publishing Company, 
New York, 1936). 


If the electron affinity which was being meas- 
ured in this experiment was essentially that of 
O., then the calculated value should be inde- 
pendent of the tube conditions. That such was 
not the case is evidenced by the fact that the 
electron affinity is seen to be markedly dependent 
on the pressure of oxygen. The average value 
obtained of 74 kcal./mole marks an upper limit 
to the electron affinity of On». 

In conclusion, it is to be noted that probably 
the main source of error is in the assumption of a 
zero reflection coefficient, but that the effect of 
a reflection coefficient different from zero would 
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be to increase still further the value obtained, 
resulting in even a wider discrepancy with the 
experimental value of Lozier’s. Also, it has been 
pointed out above that impurities, particularly 
water, can play no role since the tube was 
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adequately protected by liquid nitrogen traps. 
The formation of possible pathological molecule 
ions such as WO,- or WO;- having electron 
affinities comparable with that of atomic oxygen 
is not excluded, but is inherently unlikely. 





THE JOURNAL OF CHEMICAL PHYSICS 


VOLUME 12, NUMBER 1 


JANUARY, 1944 


The Resonance Emission of Cadmium Hydride Bands in Cadmium Photosensitized 
Reactions of Hydrocarbons* 


E. W. R. STEACIE AND D. J. LERoy 
Division of Chemistry, National Research Laboratories, Ottawa, Canada 


(Received November 13, 1943) 


An investigation has been made of the resonance excitation of cadmium hydride bands in 
mixtures of ethylene or propane with cadmium vapor illuminated with \\3261 and 2288. With 
3261 no resonance excitation with ethylene could be detected. With propane, CdH bands were 
emitted with fair intensity, thus proving that the primary process in the cadmium photo- 


sensitized reactions of the paraffins is 


Cd(5*P:)+ RH—+CdH+R. 


With 2288 no resonance emission could be detected, but the results were inconclusive. 


INTRODUCTION 


N recent years we have investigated a series of 
reactions of hydrocarbons photosensitized by 
cadmium, the two resonance lines, \3261 and 
2288, being employed.! With 3261 and the 
paraffins, the quantum yield is of the order of 0.2 
to 0.6, and the evidence points strongly to a 
C—H bond split as the initial step. On energetic 
grounds this can only occur if the primary process 
involves the formation of CdH, i.e., 


Cd(5°P1:) +RH—CdH+R, (1) 


followed by 
CdH—Cd (54S) +H. (2) 


With \3261 and olefins the quantum yield is very 
small, about 0.01, and it appears to be established 
that this is due to the quenching process involving 
the double bond. What little reaction does occur 


* Contribution No. 1174 from The National Research 
Laboratories, Ottawa, Canada. 

1See, for example, E. W. R. Steacie and R. Potvin, 
J. Chem. Phys. 7, 782 (1939). E. W. R. Steacie, D.J. 
LeRoy, and R. Potvin, J. Chem. Phys. 9, 306 (1941). 
E. W. R. Steacie and D. J. LeRoy, J. Chem. Phys, 10, 22 
(1942); ibid. 11, 164 (1943). D. J. LeRoy and E. W. R. 
Steacie, J. Chem, Phys. 10, 683 (1942), 





might be due to CdH formation by (1), or might 
involve a pseudo-thermal polymerization or some 
other mechanism. 

With \2288, reaction occurs readily with both 
olefins and paraffins. In this case the excitation 
energy is greater than the C—H bond strength, 
and it is thus possible for a C—H bond split to 
occur directly, without the intermediate forma- 
tion of cadmium hydride. 

Bender? and Olsen* have observed resonance 
excitation of the cadmium hydride bands in the 
neighborhood of \A4791, 4571, 4500, 4310, etc., 
when a mixture of cadmium vapor and hydrogen 
is illuminated with light from a cadmium-hydro- 
gen discharge, thus proving the formation of 
CdH by the process 


Cd(5°P1) +H2-CdH +H. 


It seemed of interest to carry out a similar 
investigation with hydrocarbons as a check on the 
mechanisms proposed for cadmium photosensi- 
tized reactions. 


2 P. Bender, Phys. Rev. 36, 1543 (1930). 
3L. O. Olsen, J. Chem. Phys. 6, 307 (1938), 








ca 
pre 
em 
Wi 

















CADMIUM 


EXPERIMENTAL 


The experimental method was essentially that 
employed by Bender and by Olsen. Radiation 
from a cadmium-hydrogen discharge entered a 
cell containing cadmium vapor and a hydrocarbon 
or other gas. If cadmium hydride was present in 
the cell, CdH bands were emitted as fluorescence. 

A quartz cadmium-hydrogen discharge tube 
was used, similar to that of Ellett.4 The tube was 
in the form of a spiral, and made 3} turns around 
the resonance cell. It was operated at about 2.5 
kva. It was found difficult to obtain reproducible 
results with Ellett’s arrangement in which the 
metal vapor diffused into the hydrogen stream 
from a side-tube. The arrangement was therefore 
modified so that hydrogen flowing into the dis- 
charge passed through a quartz tube heated to 
about 300°C and containing a 1 : 1 cadmium-tin 
alloy. This gave more reproducible saturation of 
the gas with cadmium vapor. The connecting 
tubing between the saturator and the discharge 
tube was wound with Nichrome ribbon and kept 
at about 350°C. 

The resonance cell was of the usual type with a 
plane window on one end and a light trap at the 
other. It was blackened except for part of the end 
window and the section opposite the discharge 
tube. It was found that a very convenient, 
adherent, flat black paint for optical purposes at 
temperatures up to 300°C can be made by 
grinding together lamp black and “‘Shellkote,” a 
high temperature lacquer. For experiments with 
3261 the resonance cell was made of Pyrex glass, 
while for \2288 fused quartz was used. The 
resonance cell and the central portion of the 
discharge tube were contained in an electric 
furnace, the heaters of which were arranged so as 
to make sure that the window of the cell was 
slightly hotter than the rest of the furnace. 

Since the purpose of the work was to investi- 
gate gases which underwent photosensitized 
decomposition, it was essential that means be 
taken to prevent effects due to the products of the 
photolysis. This was especially important in the 


case of propane, since hydrogen is one of the. 


products, and this leads readily to resonance 
emission of the CdH bands. Trial experiments 
with propane showed that sufficient decomposi- 


‘A. Ellett, J. Opt. Soc. Am. 10, 427 (1925). 
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tion occurred during the shortest feasible ex- 
posure to vitiate the results. A flow system was 
therefore used in which the contact time was 
small. The gas under investigation entered the 
system through a calibrated capillary, passed 
through a saturator containing a cadmium-tin 
alloy similar to that described above, entered the 
resonance cell near the front window, passed 
through the cell, and then on to the pump. 

The spectrograph used was a Hilger inter- 
mediate quartz spectrograph with a dispersion of 
about 75A per mm at 4500A. This was ample, 
since all that was required was the detection of 
the main cadmium hydride bands. Eastman 40 
plates were used for experiments with \3261, and 
U.V.:sensitized plates for work with \2288. 


Results with 13261 


With the flow system used it was not possible 
to obtain a reproducible supply of cadmium 
vapor when the system was evacuated, since the 
cell was open and the vapor distilled to the colder 
parts of the apparatus. It was therefore not 
practical to compare the fluorescence with the 
cell evacuated with that in the presence of a 
foreign gas. Instead a stream of nitrogen was 
used for comparison purposes, since its quenching 
efficiency is very small.® 


Nitrogen 


With nitrogen flowing through the system at 
pressures from 0.5 to 25 mm an intense green 
fluorescence was visible. With exposures of from 
10 to 60 seconds, photographs showed a very 
intense emission of the resonance line at 3261A, 
together with the main lines of the secondary 
spectrum, AA5086, 4800, 4678, 3611-14, 3466-68, 
and 3404. With long exposures the cadmium hy- 
dride bands were faintly visible at 4500 and 
4310A owing to scattered light. 


Hydrogen 


In the presence of hydrogen, in agreement with 
previous work, the whole cadmium spectrum was 
quenched, the visible green fluorescence disap- 
peared, and cadmium hydride bands at 4310, 


6 35) C. Lipson and A. C. G. Mitchell, Phys. Rev. 48, 625 
(1935). 
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4500, and 4571A were emitted with high in- 
tensity. The optimum hydrogen pressure for the 
excitation of the bands was about 0.5 mm. With 
longer exposures bands at 4791, 4835, and 4198A 
were also evident. 


Ethylene 


With ethylene at pressures from 0.5 to 7 mm, 
3261 and the secondary spectrum were strongly 
quenched as would be expected from the high 
quenching cross section of ethylene for \3261.° 
No detectable emission of the CdH bands oc- 
curred, however. In view of the very low quantum 
yield of the cadmium photosensitized reactions of 
ethylene, this is what would be expected, since a 
large formation of CdH would almost certainly be 
accompanied by a high quantum yield. The very 
small amount of reaction which does occur with 
ethylene may or may not involve the formation 
of CdH in the primary step, since the present 
method could hardly be expected to detect such 
a small concentration of cadmium hydride. 


Propane 


With propane at pressures above 1 mm ap- 
preciable quenching occurred, and CdH bands 
were emitted in fair intensity. The optimum 
pressure for the excitation of the bands was about 
7 mm. Under optimum conditions the intensity 
of the bands was of the order of } that with 
hydrogen at 0.5 mm. In view, however, of the 
relatively small quenching by propane, which 
necessitated the use of much higher pressures, it 
is not possible to make any quantitative com- 
parison of the relative efficiencies of band 
excitation by hydrogen and propane. The results 


SE. W. R. Steacie and D. J. LeRoy, J. Chem. Phys. 11, 
164 (1943). 





AND D. J. LEROY 

do indicate, however, that the excitation of the 
bands in the presence of propane is quite efficient, 
and prove definitely that at least a major portion 
of the quenching of \3261 by propane occurs by 
the mechanism 


Cd(5*P1) +CsHs—>CdH+C3H:. 


This furnishes complete confirmation of the 
assumptions which we have made previously 
regarding the mechanism of the cadmium photo- 
sensitized reactions of the paraffins. 


Results with 12288 


Experiments with \2288 were made with a 
quartz resonance cell. No attempt was made to 
filter out 43261, but instead results obtained with 
3261 and 2288 together were compared with 
those obtained with 3261 alone. Under condi- 
tions similar to those described in the preceding 
section, there was no detectable excitation of 
cadmium hydride bands by A2288 in the presence 
of either ethylene or propane. The result is rather 
inconclusive, however, since (a) \2288 was con- 
siderably weaker than \3261, (b) scattered light 
was somewhat more troublesome, (c) owing to 
the short lifetime of the 5'P; state, quenching of 
2288 by foreign gases is relatively small, and 
(d) the high absorption coefficient of cadmium 
vapor for \2288 would result in a large fraction of 
the incident radiation being absorbed near the 
wall, and wall deactivation, etc., might become 
important. It is thus not possible to decide 
whether the primary process in the cadmium 
photosensitized reactions of these substances is 


Cd(5'P:) +RH-CdH+R 
or 


Cd(5'!P1) +RH—-Cd (54S) +R+H. 











we 











THE JOURNAL OF CHEMICAL PHYSICS 








VOLUME 12, NUMBER 1 JANUARY, 1944 


An Application of Probability Considerations to the Mechanism of the Photochemical 
Polymerization of Acetylene 


BoLEestAw LuDWIK DuNICcz 
Physical Chemistry Department, Jan Kazimierz University, Lwow, Poland 


(Received April 9, 1943) 


The series of experiments described in a recent paper led to an account of the mechanism of 
the photochemical polymerization of acetylene based upon the hypothesis that this process 
should be a chain reaction operating by cycles, each of them consisting of three stages, and that 
chain termination ought to depend upon the probability of formation of an aromatic nucleus 
among the polymers resulting from each particular cycle. While nothing is known of the 
intrinsic probabilities of the reactions postulated, the proposed mechanism seemed to offer an 
opportunity for quantitative treatment by the application of probability considerations. The 
present paper shows that if the probabilities of the formation of products of the addition and 
breaking of the easily dissociated single C—C bonds in the unstable products, in the respective 
stages of reaction, are assumed equal, quantum yields and heats of reaction may be calculated, 
which are in fair agreement with data found experimentally. 





CCORDING to evidence given in the pre- 

vious paper,' in photoactivated acetylene 
in the gaseous phase there should occur polariza- 
tion of the C=C bond, and from both experi- 
mental and theoretical considerations it was 
concluded that the addition of a normal acetylene 
molecule to an activated one, following the pre- 
vious association of these molecules, should give 
an activated cyclobutadiene molecule. No evi- 
dence was found suggesting the occurrence of any 
direct dimerization of acetylene molecules to 
vinylacetylene. 

Similarly, the first stage of any other cycle of 
this reaction would consist of the association of a 
normal acetylene molecule to the activated 
C=C bond regenerated in the growing polymer, 
followed by the formation of the activated cyclo- 
butadiene ring. 

Cyclobutadiene, or its derivatives, has not, 
so far, been obtained experimentally in any 
traceable amounts, this being usually attributed 
to the great strain to be expected, for stereo- 
chemical reasons, in the cyclobutadiene ring. 
But Wheland,? by use of the molecular orbital 
treatment, deduces that the cyclobutadiene 
ring in an activated state (i.c., lowest triplet 


state) ought to manifest properties analogous | 


to those of the benzene ring in the ground state. 
This also justifies our expecting the reactions 


1B. L. Dunicz, J. Am. Chem. Soc. 63, 2461 (1941). 
* Wheland, Proc. Roy. Soc. (London) A164, 397 (1938). 


of activated cyclobutadiene to be similar to the 
reactions of benzene in the ground state. The 
addition, therefore, of an acetylene molecule to 
the activated cyclobutadiene molecule* should 
not cause any primary dissociation of the cyclo- 
butadiene ring,** just as the benzene ring does 
not dissociate immediately in any reaction of the 
benzene molecule, but only as the result of some 
previous addition reaction to the ring. There 
would then exist two ways of acetylene reacting 
with a cyclobutadiene molecule: 


A: By association of acetylene by its C—H bond 
to cyclobutadiene: 

1. Assuming for cyclobutadiene a structure 

like the Armstrong-Baeyer structure for benzene: 


CH—CH H CH—CH, 

)» a m | = | | 

ran i 

CH—CH ‘=CH CH—CH—C=CH 


acetylene-1, 
cyclobutene-2. 


2. Assuming the possible isomeric cyclobutine: 
CH:—C H CH:—CH 

| it | 

CH,—C C=CH CH,—C—C=CH 


acetylene-1, 
cyclobutene-1. 


| 


*In the subsequent discussion, since cyclobutadiene is 
always assumed to be formed in its activated state, no 
further description of its state will be made. 

** Otherwise one would expect aliphatic trimers to be 
formed in the course of the reaction. But this is contrary 
to the experimental evidence. 
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B: By association of acetylene by its C=C bond 
to cyclobutadiene: 

1. Assuming for cyclobutadiene a structure 
like the Armstrong-Baeyer structure for benzene: 


TaBLeE I. Likely and unlikely isomerizations of the un- 
stable acetylene trimers J, JJ, III, involving breaking of 
a single C—C bond followed by intramolecular migration 
of not more than one hydrogen atom. The bond under 
consideration is dotted. 


Unlikely products. 
Their structures 
should involve 
some further 
rearrangements 
into aliphatic 
compounds, that, 





Unstable trimer of 
acetylene undergo- 
ing isomerization 





however, have H-atom H-atom 

not been found migra- migra- : 

experimentally. tion < _ tion Likely products 
I 


CH—CH CH: 68 


| | | 
CH—C——CH 


B 
CHC ois ary 


oo 


with two stereo- 
chemically dif- 
ferent types of 
C—C bond 


¥ 
i wae « 
CH—CH—CH 
a 6 a 
Bi i B 
aie “oH 
| 
CH—CH—CH 
ay 6 a 
II 


with three stere- 
ochemically dif- 
ferent types of 
C—C bon 


Bi Y 
oe 
| {| 
(m—C~CH 
a 6 a 


1 + £ 
ere 
H,—C—CH 

a 6 a 

1 +. 2 
CH2,—C—CH 


ay 


os 


All five single 
C—C bonds in 
I are likely to 
break giving in 
one case aro- 
matic ring in 
other four the 
acetylene deriv- 
ative. 


we as 
| 
CH—CH=CH 
benzene 
a? ws CH 
CH—CH—C 


acetylene-1, 
cyclobutene-2 


Only two single 
C—C bonds of 
the a, B, , 4, 
ring in JI are 
likely to break 
giving the acet- 


ylene _ deriva- 
tive. 

nit % ii 
CH.—C——C 


acetylene-1, 
cyclobutene-1 
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TABLE I[.—Continued 











Unlikely products. 
Their structures 
should involve 

some further 
rearrangements 

into aliphatic 
compounds, that, 


Unstable trimer of 
acetylene undergo- 
ing isomerization 


however, have H-atom H-atom 
not been found migra- migré a- 
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with five stereo- 
chemically dif- 
ferent types of 





Only three sin- 
gle C—C bonds 
of the a, B, vy, 4, 




















C—C bon ring in JJI are 
likely to break 
giving in one 
case aromatic 
ring, in other 
two the acety- 
lene derivative. 

Bi ¥ 8 

CH.—CH — Bid saan —CH 

| : — | 

CH =C——CH CH=CH—CH 

a1 6 a benzene 

Bi ¥ B 

i. 7 _ ya oo. B5 lle . 
| -— | = III 
Cn,.—C—On CH=C---CH CH—CH—C 
a1 6 a acetylene-1, 
cyclobutene-2 

Bi B 

CH, —CH- -CH ay ‘a a 

|| —> 
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a 6 a acetylene-1, 
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8 ¥ B 
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2. Assuming the possible isomeric cyclobutine: 





CH 
+ || os 
CH 


CH.—C—CH 
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CH.— 
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3. The formation in either of the above reac- 
tions (B1 and 2) of the third trimer should also 
be possible: 





1 2 
CH.—CH—CH 
III. | . be | 
CH =C——CH. 





4 


Owing to the undoubtedly unstable strained 
structure of the acetylene trimers J, JJ, III, it 
would be difficult to imagine that the process 
of adding acetylene molecules could continue 
beyond these products, without previous iso- 
merization into some more stable forms occuring. 
In Table I are shown the isomerizations which 
might be involved. They are based on the as- 
sumption of only one single C—C bond breaking, 
this being followed by intramolecular migration 
of not more than one hydrogen atom. The 
products considered as likely would then be: 


acetylene-1, cyclobutene-1, 
acetylene-1, cyclobutene-2, 
and benzene.t 


In Table I the respective products of the iso- 
merization reactions considered have been classi- 
fied as likely and unlikely. Reactions are as- 
sumed to be likely when the breaking of a single 
C—C bond gives rise to an aromatic ring, or to 
an acetylene-cyclobutene structure with regener- 
ation of an activated C=C linkage capable of 
starting a new cycle of reactions by forming a 
cyclobutadiene ring with a new acetylene mole- 
cule. Reactions are considered unlikely if their 
products cannot be regarded as _ sufficiently 
strainless not to undergo further rearrangements 
to attain greater stability—such products never 
having been observed experimentally. 

As the type A of reaction does not predict the 
possibility of benzene formation, in contradic- 
tion to the experimental findings of various in- 
vestigators mentioned elsewhere,! one is tempted 
to conclude that the reaction which occurs is of 


the B type. This could also be deduced by con- 


sidering that the triple bond is, because of its 


+ This stands in agreement with experimental observa- 
tions of the formation of benzene and naphthalene in the 
photochemical polymerization of acetylene pure and with 
added benzene, respectively, in the gaseous phase (see 
reference 1). 





OF C:H: 39 


large accumulation of electrons, more liable than 
a homopolar, single, C—H bond, to polarization 
by induction. This makes possible the conse- 
quent association of acetylene, by its C=C bond, 
to the cyclobutadiene molecule. If then we as- 
sume the B type of reaction as being the most 
likely, the simplest further assumption we can 
make is that the various fused-ring (bicyclic) 
structures which can be formed in this reaction, 
have equal chances of formation. Similarly in 
these fused-ring structures themselves, we as- 
sume equal chances for the breaking of the single 
C—C bonds capable of giving rise to those 
products of isomerization classified as ‘‘likely.” 
The only satisfactory alternative to this would 
be to attempt to evaluate the intrinsic probabili- 
ties of the various reactions involved, by the 
kinetic treatment. But this, if possible at all, 
would be very complex. 

With regard to the three unstable trimers: 
I, II, IIT, of acetylene, to be expected from the 
first cycle of reaction, there have been ascribed 
therefore equal chances for the dissociation of 
all five single C—C bonds, marked 1-5 with the 
probability 4 for each, for the case of trimer J; 
equal chances for the dissociation of two bonds, 
marked 2 and 4, with the probability 4 for each, 
making 1 for the whole system (as breaking either 
of these bonds would give the same result) for 
the case of trimer JJ; and equal chances for the 
dissociation of three bonds, marked 2, 4, and 5, 
with the probability 4 for each, for the case of 
trimer JJI. Similar conditions have been as- 
sumed to apply to the successive cycles of the 
photochemical change of acetylene and also to 
the reactions of the mono- and di-substituted 
acetylenes, and on such considerations there have 
been worked out Tables III-VIII, containing 
totals of the probabilities (i.e., total yields) of the 
reactions occuring in particular cycles of the 
respective polymerizations. A better insight into 
Tables III-VIII is given by Table II, where is 
shown the probable effect of substituents on the 
breaking of the single C—C bonds in the un- 
stable fused-ring structures produced in the 
second stages of the cycles. Breaking of the 
single C—C bond is considered improbable if the 
substituent replaces a hydrogen atom capable of 
migrating so as to complete isomerization into 
likely products, i.e., either a deactivated polymer 


(containing an aromatic ring) or a polymer cap- 
able of further reaction by virtue of possessing 
an activated C=C bond regenerated in it. 

From the Tables III-VIII it was then possible 
to derive the following means of calculating 
quantum yields and heats of reaction for several 
acetylenes: 


For the case of acetylene, one might illustrate | 


schematically the whole course of its polymeriza- 
tion, as shown in Tables III-V, by depicting the 
course of the reaction assumed to take place in 
the first three of its cycles. In the formulae used 
in these tables the single, double, and triple 
intercarbon bonds are presented by means of 
single strokes, and double and triple strokes 
joined at the ends, respectively, carbon atoms 
being understood to be placed at their ends. The 
hexagon and the small bisected square denote 


TABLE II. Effect of substituents on the breaking of the 
single C—C bonds in the unstable fused-ring structures 
produced in the second stages of the cycles. (The single 
C—C bonds capable of breaking are shown by dotted lines. 
Substituents which restrict bond breaking are represented 
by R and a number. Those unlikely to influence bond 
breaking are represented by R*. Not more than one sub- 
stituent is assumed to be possible on any carbon atom. 
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TABLE II.—Continued 
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benzene and cyclobutene rings, respectively. For 
simplicity the bonds between hydrogen and 
carbon atoms are not shown; their number and 
positions can be traced from the number of repre- 
sented bonds meeting together, or from compar- 
ing the respective bond arrangements in these 
formulae with the formulae in Table II. In the 
interstices Ao, Ai, A2, there are shown the first 
stages of reaction in the respective cycles: forma- 
tion of an activated cyclobutadiene ring from two 
C=C linkages, one of them being in an activated 
state. In the interstices Bo, Bi, Bo, there are 
shown all the possibilities, allowed by the postu- 
lated mechanism of reaction, for the formation of 
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the unstable products of the reaction between 
cyclobutadiene and the next acetylene molecule. 
These possible structures are derived by con- 
sidering the combination of the substituents with 
the bonds in the cyclobutadiene/cyclobutine/ 
ring, so as to give all possible structural isomers. 
Differences between the substituents are ignored, 
the only condition required being that the asso- 
ciated acetylene molecule has one link in the 
C=C bond opened in the addition reaction. The 
yield (probability of formation) of any of such 
products is taken as the reciprocal of the number 
of all possible structural isomers, permitted by 
the mechanism. In the interstices Co, Ci, Cs, 
there are two columns showing in what propor- 
tion the unstable products formed in the second 
cycle would undergo isomerization into two 
classes of stable products: (a) those deactivated 


TABLE III. The first cycle of the photochemical poly- 
merization of acetylene. Ao—Ist stage of the cycle: forma- 
tion of an activated cyclobutadiene ring, non-, dimethyl-, 
or tetramethyl-substituted, from two acetylene, methyl- 
acetylene, or dimethylacetylene molecules, one of them 
being activated, respectively. (In case of methylacetylene 
two dimethylcyclobutadienes are possible; a yield of 1/2 
for each is assumed.) Bo—2nd stage of the cycle: bonding of 
acetylene, methylacetylene, or dimethylacetylene molecule 
to the respective activated cyclobutadiene ring, formed in 
the first stage. The postulated mechanism of reaction would 
allow for primary formation of stereochemically different, 
unstable, fused-ring structures: three in case of acetylene 
with assumed yield 1/3 for each; twelve in case of methy]l- 
acetylene with assumed yield 1/12 for each; and one in 
case of dimethylacetylene. Co—3rd stage of the cycle: the 
unstable fused-ring structures formed in the second stage 
undergo rearrangement into either an aromatic nucleus 
(first column) causing deaciivation of the product, or an 
acetylene-cyclobutene structure (second column) with 
activated C=C bond regenerated, giving start to the 
second cycle of reaction. The derivation of the respective 
probabilities of rearrangement, shown by fractional values, 
is to be found in the text. The general result for the whole 
cycle is that from three moles of acetylene, methyl-, or 
dimethyl-acetylene, respectively, one mole being activated, 
there is formed: in the yield ap—the deactivated trimer, 
i.e., One containing an aromatic nucleus; in the yield bo>— 
the reactive trimer containing a regenerated activated 
C=C bond, giving start to the second cycle of reaction. 
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through formation in them of an aromatic ring,* 
and (b) those undergoing further reactions 


TABLE IV. The second cycle of the photochemical poly- 
merization of acetylene. Ai—Jst stage of the cycle: forma- 
tion of an activated cyclobutadiene ring, mono-, or tri- 
substituted, through bonding of acetylene or methyl- 
acetylene molecule to the activated C==C bond of the 
respective reactive trimer of acetylene or methylacetylene, 
resulting from the first cycle. Bij—2nd stage of the cycle: 
successive bonding of acetylene or methylacetylene mole- 
cule to the respective activated cyclobutadiene ring, formed 
in the first stage. The postulated mechanism of reaction 
would allow for primary formation, in each respective case, 
of six stereochemically different, unstable, fused-ring struc- 
tures with yield 1/6 assumed for each (difference between 
methyl- and dimethylcyclobutenyl-substituents in case of 
methylacetylene reaction being neglected). C:;—J3rd stage 
of the cycle: the unstable fused-ring structures formed in 
the second stage undergo rearrangement into either an 
aromatic nucleus (first column) causing deactivation of the 
product, or an acetylene-cyclobutene structure (second 
column) with activated C=C bond regenerated, giving 
start to the third cycle of reaction. The derivation of the 
respective pobabilities of rearrangement, shown by frac- 
tional values, is to be found in the text. The general result 
for the whole cycle is that from one mole of reactive trimer 
(of acetylene or methylacetylene), resulting from the first 
cycle, together with two moles of acetylene or methyl- 
acetylene, there is formed: in the yield a—the deactivated 
pentamer, i.e., one containing an aromatic nucleus; in the 
yield b—the reactive pentamer containing a regenerated 
activated C=C bond, giving start to the third cycle of 
reaction, in which the yields of the deactivated and re- 
active polymers will again be a and 3, respectively. 
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* This is in accord with the observation that the con- 
centration of benzene steadily increases during the reaction 
while that of other volatile products attains some sta- 
tionary equilibrium concentration, whose value depends 
upon the kinetics of the reaction. If, on the other hand, 
deactivation of the growing polymer were not exclusively 
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TABLE V. The third cycle of the photochemical poly- 
merization of acetylene. (The case of a reactive C=C bond 
between two cyclobutene groups; the yields a and b of 
acetylene polymers, deactivated and undergoing the 
further reaction, respectively, should be the same as in 
case of cyclobutene groups at one end of the reactive 
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through regeneration in them of the activated 
C=C linkage, these products being multiplied 
by the probability of formation of the unstable 
product of the second cycle. Their sums a» and 
bo or a and b, shown at the bottom of the tables, 
are the probabilities of formation in the cycle 
of the two respective classes of products. 

In the case of methyl substituted acetylenes it 
is necessary to take into account the possibility 
that these compounds when activated, might 
undergo isomerization giving allenes, probably 
through intermediate formation of diradicals, e.g., 
—CH.,—CH=CH-—-, —CH2.—CH=C(CHs)-. 
In the photochemical reactions of these com- 
pounds, therefore, twofold reactions might be 
expected to occur: (a) A reaction similar to that 
of acetylene, as shown in Tables VI-VIII (the 
symbolic formulae in these tables have the same 
significance as in previous tables with addition 


associated with the formation of an aromatic nucleus 
(according to this derivation), one would expect to find 
also other volatile aromatic products, such as styrene or 
naphthalene, which, however, have not been experimen- 
tally observed. 





TABLE VI. The first cycle of the photochemical poly- 
merization of methylacetylene reacting in a similar way to 
acetylene. Some difference in further reacting of trimers 
with reactive C=C bond having substituents at one and 
both ends may be neglected on account of the very small 


probability of the former). 
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of black circles (e) to represent methyl groups, 
and also in Table VII dimethylcyclobuteny]l), 
with probabilities of formation of the two classes 
of products, deactivated and undergoing further 
reaction, different from those of the change of the 
acetylene, due to substituents. (b) Formation 
from the methylacetylene of an intermediate 
diradical, which then becomes deactivated by 
reacting with a second molecule of acetylene 
giving a hydroaromatic cyclopentadiene com- 
pound. This is shown in Table IX. 

In the case of deuterium substituted acetyl- 
enes one would expect reduced mobility of the 
deuterium atom in the molecule as compared 
with the lighter hydrogen atom; monodeutero- 
acetylene should therefore show a yield inter- 
mediate between that for acetylene and that for 
methylacetylene undergoing type (a) reaction 
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(i.e., similar to that of acetylene), and dideutero- 
acetylene ought to show a still greater tendency 
to termination of reaction chains through the 
formation of aromatic rings. 

Assume now the reaction to be started by one 
mole of activated acetylene molecules. As shown 
in Tables III, VI, and VIII, 3 moles of acetylene 
would be involved in the first cycle. Let ao be 
then the probability of formation of the aromatic 
nucleus from the three C=C linkages of the 
respective acetylenes. It follows that the number 
of acetylene molecules used up in the formation 
of the deactivated polymer will be 3a9; at the 
same time the heat evolved by the formation of 
this product would be approximately 149d» kcal. 
in this cycle.* In the formation of products 
capable of further reaction, the number of 
acetylene molecules used up would be 3bo 
(where b)=1-—ao); and the heat evolved will 
be approximately 88d» kcal.** The 3b acetylene 


TABLE VII. The second cycle of the photochemical poly- 
merization of methylacetylene reacting in a similar way to 
acetylene. (Black circles denote methyl groups and also the 
dimethyl-cyclobutenyl groups which originated from the 
first cycle.) 
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* The heat of formation of the aromatic nucleus from 
three C=C bonds is taken as the difference between the 
heat of combustion of three moles of acetylene (3X 312.4 
= 937.2 kcal./mole) and one mole of benzene vapor (788 
kcal./mole), i.e., 149.2 kcal. Seé Landolt-Bérnstein, Physi- 
kalisch-chemische Tabellen (Berlin, 1923), series 5, Vol. 2, 
pp. 1588-9. 

** Of the three C=C bonds, with bond energy 123 
kcal./mole each, reacting to form an acetylene-cyclo- 
butene, two are being transformed into one C=C bond, 
with bond energy 100 kcal./mole, and four C—C bonds, 
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molecules used up in the formation of the reac- 
tive products, would combine in the second 
cycles, as shown in Tables IV and VII, with an 
additional 2b) acetylene molecules. Hence for the 
polymer resulting from the second cycle the 
total number of acetylene molecules used up 
would be (3+2)o, of which (3+2)boXa would 
produce the new deactivated polymer (pen- 
tamer), the heat of reaction being approximately 
(149+88)boXa kcal.*** The number of moles 
of acetylene in the reactive products, resulting 
from the second cycle, will be (3+2)boxéb 
(where b=1—a). These products would react 
further in the third cycle. It may be observed 


TABLE VIII. The photochemical polymerization of di- 
methylacetylene reacting in a similar way to acetylene. 
(From three moles—one activated—of dimethylacetylene 
molecules involved, all should yield hexamethylbenzene, 
this being deactivated.) 


| det fe] xe 
| net rH] | - 


@o=t. De= 0 















































TABLE IX. The course of the photochemical polymeriza- 
tion of methyl- and dimethyl-acetylene, reacting through 
the formation of diradicals:s —CH2—CH=CH—, and 
—CH:—CH=C(CH;)—, respectively, upon activation. 
The products would be hydroaromatic dimers: methyl- and 
trimethyl-cyclopentadiene, respectively, these being de- 
activated. 
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with bond energy 58.6 kcal./mole each, the third C=C 
being regenerated. Thus the heat of formation of the 
acetylene-cyclobutene system from acetylene would be 
the difference between the sum of bond énergies of 1 C=C 
and 4 C—C bonds and the bond energies of 2 C=C bonds, 
giving 88.4 kcal./mole of the product. See Pauling, The 
Nature of the Chemical Bond (Cornell University Press, 
1942), second edition, pp. 53 and 131. 

*** The product comprises the benzene-cyclobutene 
system formed from five C=C bonds, three of them giving 
the benzene ring with evolution of 149 kcal./mole of 
benzene. The other two form the cyclobutene ring and the 
bond between benzene and cyclobutene rings, the heat of 
formation evolved being 88 kcal./mole. 
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TABLE X. Progress of formation of deactivated polymers in the photochemical polymerization of acetylene, 
started by a mole of activated acetylene molecules. 








Number of acetylene 


Heat of reaction evolved in each 
cycle on formation of a deactivated 





Num- molecules in the Number of acetylene molecules polymer, in kcal./mole 
ber of polymer resulting producing a deactivated polymer of activated acetylene 
cycle from the cycle in each cycle molecules involved 
1 3 39 149a9 
2 5 tt boXa Hetenties 
3 7 3+2.2boXbXa [149+88.2 JboxbXa 
4 9 [3+2.3 looXb?Xa [149+88.3 ooXb?Xa 
n—1 3+2(n—2) PE aS coy io Hie eote— 1h pace 
n 3+2(n—1) 3+2(n—1) ]boXb"2xXa 149+88(n—1) ]boXb"?*Xa 








from the tables that each successive cycle would 
result in two additional acetylene molecules being 
brought into reaction in stoichiometric relation, 
as outlined in the Table X. Hence in the forma- 
tion of the polymer resulting from the third cycle 
the total number of acetylene molecules in- 
volved would be (3+2.2)boXb. The number of 
moles of acetylene in the deactivated polymer 
produced in the third cycle will be (3+2.2)do 
Xb Xa, and the evolved heat on its formation be 
approximately (149+88.2)bo Xb Xa kcal. 

For the cycles following the second one, the 
probability of formation of the deactivated 
polymer in each cycle would be always a(as in 
the second cycle), and the number of acetylene 
molecules forming a reactive polymer remaining 
after each cycle would be then } times the number 
of molecules contained in the reactive polymer 
resulting from the previous cycle, after the 
addition to it of 2 acetylene molecules in stoichi- 
ometric relation. Therefore, for the polymer 
resulting from the fourth cycle the number of 
moles used up would be (3+2.3)b9X0?, and the 
number used up for the formation of the de- 
activated polymer would be (3+2.3)bo Xb? Xa, 
the heat evolved being (149+88.3)b) xb? Xa 
kcal. In the same way it may be shown that 


+ From the second cycle onwards the heat evolved in the 
formation of the deactivated polymer would be made up 
as follows: 

(a) The heat of formation of the cyclobutene rings con- 
tained in it multiplied by the probability of formation 
of the polymer. The number of cyclobutene rings will 
be one less than the number of cycles, hence the heat 
evolved will be 88(m—1) kcal. 

(b) The heat of formation of one aromatic nucleus termi- 
nating the chain of reaction, i.e., 149 kcal./mole. This 
is also multiplied by the probability of formation of 
the polymer. 

The total heat of formation will be the sum of these two 


terms. 





for the deactivated polymer of the mth cycle 
the number of molecules used up would be 
(3+2(n—1))bo Xb" Xa and the heat evolved on 
its formation (149+88(n—1))bo Xb" Xa kcal. 

Hence the total number of acetylene moles 
involved in the polymerization initiated by one 
mole of activated acetylene molecules, with 
cycles running into infinity, would be: 


g=3ao+ r [3+2(n—1) ]boxXb"-*Xa 


n=2 


= 3a9+3aXbo > b"+2a 


m=0 


Xbo > (m+1)b" (putting n—2=m) 








m=0 
=3a9t+3aXboX +2aX bo 
x (by the binomial theorem) 
(1—5)? 
2bo 
= 3a) +3bo+— 
a 


=3+2b)/a. 


This last quantity represents the theoretical 
quantum yield ¢ of the reaction. 

Correspondingly, the heat evolved in the 
reaction would be: 


149a.+ > [149+88(n—1) ]bo Xb"? 
n=2 
Xa=149+88b)/a kcal. 


To obtain the heat change per acetylene mole 
polymerized, this quantity must be divided by 
the quantum yield g; thus the theoretical heat 
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TABLE XI. Data for the photochemical and radiochemical polymerization of acetylene, 
methyl-, dimethyl-, deutero-, and dideutero-acetylene. 











Compound 


Q kcal./mole of monomer 





P calc. P exp. M/2N exp. calc. exp. 

CH=CH P=9.8 9.2+1.5 (6) 10-13 (5, 9, 10) Q2>45.7 47 (10)* 
CH;—C=CH 2 <€P<37 3.3 (7) 4.15 (2) 47.1>Q>42.5** —_ 
CH=CD 5.7<P<9.8 —_ —_ 47.1>0>45.7 ~ 
CH;—C=C—CH; 2 €P<3 — 2.9 (2) 49.3>Q>42.5** — 
CD=CD 3 <P<5.7 5 (3) om 49.3>0>45.7 _ 








* For radiochemical polymerization. 


** Heat of formation of cyclopentadiene from acetylene and methylacetylene, on the basis of available data. See Landolt-Bérnstein, Physicalish- 
chemische Tabellen (Berlin), series 5, Vol. 2, pp. 1588-9 (1923), and Vol. 3, p. 2896 (1936). 


of reaction would be: 
Q=(149+880)/a)/¢ kcal. 


Since in the case of methyl and dimethyl 
acetylenes the chain of reaction might be stopped 
if these were to dimerize partially into methyl- 
and trimethyl-cyclopentadiene, respectively, the 
heat of polymerization of these compounds per 
mole acetylene polymerized should lie between 
the values calculated for their reaction in a 
similar way to acetylene, and one-half of the 
heat of formation of the cyclopentadiene ring 
from one C=C and one C=C—C linkage, i.e., 
approximately 42.5 kcal./mole of acetylene 
polymerized.* 

The final values of the quantum yield and 
heat of reaction for various acetylenes considered, 
as calculated by this method, are shown in 
Table XI, where they are compared with the 
available experimental data. In the case of the 
heats of reaction any additional heat which 
might be evolved in the secondary combination 
of unsaturated polymers could not be taken into 
account here, although this quantity might be 
appreciable and would make the heat of reaction 
more exothermic. As may be seen, there is good 


* The heat of formation of the cyclobutadiene ring has 
been taken as the difference between the sum of heats of 
combustion of acetylene and methylacetylene (312.4 kcal./ 
mole of acetylene+473.6 kcal./mole of methylacetylene 
= 786.0 kcal.) and the heat of combustion of liquid cyclo- 
pentadiene (701.0 kcal./mole, this being the only available 
value). See Landolt-Bérnstein, Physicalisch-chemische 
Tabellen (Berlin), series 5, Vol. 2, pp. 1588-9 (1923), and 
Vol. 3, p. 2896 (1936). 





agreement between the experimental and the 
calculated values for all the compounds con- 
sidered. The data suggest also that both the 
radiochemical and photochemical reactions fol- 
low a similar mechanism. 

The somewhat larger yields in the radiochem- 
ical reactions of acetylenes under the influence 
of both alpha- and beta-rays** might possibly 
be due to the fact that ionization of acetylene 
molecules might be accompanied by polarization 
of some number of normal molecules. Moreover, 
because of the enormous energy released by the 
disintegration of radon,’ one ought to expect a 
larger yield of aromatic products®*-” in the 
alpha-ray polymerization of acetylene than in the 
photochemical reaction, since it is well known 
that the effect on acetylene of a rise in tempera- 
ture is to favor the formation of aromatic hydro- 
carbons. This should also provide an explanation 
for the observed temperature independence of 
the yield of the radiochemical reaction** within 
the range 0°-100°C. 


3 Heisig, J. Am. Chem. Soc. 53, 3245 (1931). 
a : o Bardwell, and Perry, J. Am. Chem. Soc. 48, 1556 
“a Mund and Jungers, Bull. Soc. Chim. Belg. 40, 158 
31). 
6 Mund, Velghe, Devos, and Vanpee, Bull. Soc. Chim. 
Belg. 48, 269 (1939). 
7 Rutherford, Chadwick, and Ellis, Radiations from 
Radioactive Substances (Cambridge, 1930), p. 162/3. 
8 Rosenblum, Bull. Soc. Chim. Belg. 46, 503 (1937). 
® Jungers and Taylor, J. Chem. Phys. 3, 338 (1935). 
10 Lind and Livingston, J. Am. Chem. Soc. 54, 94 (1932). 
( ul oi and Livingston, J. Am. Chem. Soc. 55, 1036 
1933). 
12 Lind and Schiflett, J. Am. Chem. Soc. 59, 411 (1937). 


